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Summary 


The  project,  “Near-Field  Propagation  of  Sub-Nanosecond  Electric  Pulses  into  Amorphous 
Masses”,  deals  with  the  propagation  of  subnanosecond  pulsed  radiation  in  the  near  field  region 
and  its  use  in  probing  inhomogeneities  in  dielctrics.  It  also  explores  the  possibility  of  changing 
the  physical  structure  of  lossy  dielectrics  by  using  very  intense  pulses.  Our  research  effort  has 
been  focused  on  five  research  areas: 

•  Design  and  construction  of  a  compact  subnanosecond  pulsed  power  generator 

•  Focusing  of  subnanosecond  pulsed  electric  fields  in  the  near  field  using  prolate  - 
spheroid  wideband  antennas 

•  Increase  of  the  spatial  resolution  (reduction  of  the  focal  volume)  using  a  dielectric  lens  in 
addition  to  the  reflective  antenna 

•  Target  detection  with  focused  subnanosecond  radiations 

•  Modification  of  structure  of  biological  targets  with  intense  subnanosecond  pulses 

A  Subnanosecond  High  Voltage  Pujse  Generator 

A  subnanosecond  high  voltage  pulse  generator  was  constructed,  which  is  able  to  deliver  150 
ps  (FWHM)  pulses  with  100  ps  risetime  (10%-90%).  The  maximum  voltage  is  17  kV  and  the 
repetition  rate  is  20  Hz.  This  pulse  generator  uses  a  four-stage  Marx-Bank  to  charge  a  peaking 
capacitor,  which  then  discharges  through  a  peaking  switch  into  a  transmission  line  [1-10].  As  a 
result,  a  step  voltage  with  a  fast  risetime  followed  by  an  exponential  decay  is  generated.  The 
step  voltage  is  then  cut  off  by  a  tail-cut  switch  to  achieve  the  desired  pulse  durations.  A 
photograph  of  the  pulse  generator  is  shown  in  Fig.  la. 

The  fast  rise  and  fall  time  of  the  pulses  generated  with  this  system  is  due  to  the  fast  voltage 
rise  of  the  pulse  generated  by  the  Marx  generator  (6.6  kV/ns).  Because  of  the  high  dV/dt,  the 
voltage  across  the  peaking  switch  (40  kV)  reaches  values  of  almost  eight  times  the  static 
breakdown  voltage.  The  high  overvoltage  is  the  reason  for  the  extremely  fast  breakdown  in  the 
temporal  range  of  100  ps.  The  tail-cut  switch  which  is  integrated  into  the  coaxial  cable  in  front  of 
the  load  is  similarly  over-volted,  leading  to  a  fast  cut-off  of  the  pulse  on  a  time  scale  of  again 
100  ps.  The  temporal  development  of  the  pulse  is  shown  in  Fig.  1b. 

Unlike  most  of  the  subanosecond  switches,  that  are  operated  in  gases  (in  hydrogen  for 
example)  at  a  pressure  several  times  higher  than  atmospheric  pressure,  the  spark  gap  switches 
in  this  pulse  generator  are  operated  in  atmospheric  pressure  air.  This  allows  for  the  construction 
of  a  simple  and  compact  pulsed  power  generator  (Fig.  la). 

In  addition  to  this  pulse  generator  a  commercially  available  pulse  generator  which  provided 
200  ps  pulses  at  an  amplitude  of  up  to  5  kV  was  used  for  the  studies  [FPG5-10PM,  FID]. 
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Fig.1.  a)  Photograph  of  the  four-stage  Marx  generator  including  peaking  switch.  For  size-comparison  a 
calculator  is  shown;  b)  Temporal  development  of  a  16  kV,  150  ps  (FWHM)  pulse  across  a  50  Q  load. 


Focusing,  subnanosecond  pulses  in  the  near  field 

Focusing  electromagnetic  waves  into  soft  matters  were  conducted  using  annular  antenna  arrays 
[11-17],  Narrow  band  signals  were  emitted  from  the  array  antennas  and  created  a  constructive 
interference  at  the  desired  focal  point.  Ultrawideband,  pulsed  signals  were  studied  using 
delayed  beam  forming  [18]  and  time-reversal  technique  [19-20].  These  two  techniques  were 
limited  to  low  power  signals,  which  allow  for  high-accuracy  temporal  coordination  among 
antennas.  High  power,  narrowband  reflector  antennas  for  focusing  in  the  near  field  were  studied 
in  [21-23].  The  focusing  cavity  was  even  filled  with  low-loss  dielectrics  and  metamaterial  in  order 
to  achieve  a  better  transmission  from  the  antenna  to  the  matter  [24-27],  The  high  power, 
subnanosecond  pulsed  radiation  was  focused  by  means  of  a  prolate-spheroidal  impulse 
radiating  antenna  (IRA),  which  was  developed  by  Carl  Baum  at  the  University  of  New  Mexico 
[28-31],  and  built  by  Farr  Research  (Albuquerque,  NM)  (Fig. 2).  In  this  antenna,  spherical  waves 
emanate  from  the  first  focal  point  of  the  prolate-spheroidal  reflector  and  propagate  along  two 
pairs  of  parallel,  balanced  conical  transmission  lines.  The  electromagnetic  waves  related  to  the 
subnanosecond  pulse,  that  are  scattered  on  the  reflector  surface  and  converge  to  the  second 
focal  point,  are  called  “impulse”.  Another  part  of  the  spherical  waves  which  is  propagating 
opposite  to  the  impulse  is  called  “prepulse”  as  it  arrives  at  the  second  focus  earlier  than  the 
impulse.  The  prepulse  will  not  be  focused  and  its  amplitude  decreases  with  distance. 

The  amplitudes  of  both  prepulse  and  impulse  are  plotted  in  lateral  and  (longitudinal  (axial) 
direction  in  Fig.  3  [32],  The  power  density  of  the  impulse,  which  is  proportional  to  the  square  of 
the  electric  field,  is  focused  in  both  lateral  and  axial  directions.  The  distribution  of  power  density 
has  a  FWHM  width  of  32  cm  in  the  axial  direction  and  a  FWHM  width  of  10  cm  in  the  lateral 
direction  (at  z=16  cm). 
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Fig. 2.  Antenna  structure,  (a)  Side  view.  A  differential  signal  is  fed  to  the  transmission  line  at  the  first  focal 
point  (-z0).  (b)  Front  view  <t>0  is  the  angle  between  the  feed-arm  and  y-axis.  (c)  The  actual  antenna. 


(a)  (b) 


Fig. 3.  (a)  Measured  electric  field  (x-directed)  distribution  along  the  axis  (z-axis).  The  prepulse,  the  impulse, 
and  the  square  of  impulse  are  shown.  The  curve  with  down-triangle  symbols  is  the  modeled  impulse 
distribution  along  the  z-axis.  (b)  The  measured  electric  field  (x-directed)  distribution  in  the  lateral  direction 
(along  the  y-axis)  at  the  focal  plane.  The  prepulse,  the  impulse,  and  the  square  of  impulse  are  shown. 
The  curve  with  down-triangle  symbols  is  the  modeled  impulse  distribution  in  y-axis. 

The  electric  field  distribution  of  this  antenna  was  modeled  by  using  CST  Microwave  Studio,  a  3- 
D,  finite  integral  time  domain  (FITD),  electromagnetic  simulation  software.  The  modeled  results, 
for  emission  into  free  space  agree  very  well  with  the  experimental  results.  This  holds  especially 
near  the  focal  point  for  both  the  axial  and  the  lateral  electric  field  distributions.  Closer  to  the 
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reflector  edge  we  observe  a  slight  discrepancy  between  the  modeled  field  and  measured  field, 
which  could  be  due  to  the  fact  that  the  model  only  allows  us  to  compute  the  x-directed  field,  but 
the  sensor  may  pick  up  electric  field  components  in  y  and  z  direction. 

The  spatial  distribution  of  the  electric  fields  in  the  near  field  of  the  focusing  antenna  were  also 
modeled  for  dielectric  media,  other  than  in  air.  Because  the  radiation  from  the  reflector  antenna 
is  composed  of  spherical  waves  converging  to  the  second  focus,  a  hemispherical  shape  was 
chosen  for  the  dielectric  with  its  center  at  the  second  focal  point.  In  this  case  the  phases  of 
incident  waves  at  the  periphery  of  the  hemisphere  are  equal.  In  the  simulation,  the  dielectric 
hemisphere  has  a  diameter  of  12  cm  and  dielectric  constant  (er)  of  12.  Different  conductivity  (a) 
values,  0,  0.3  S/m  and  0.5  S/m  were  assigned  to  the  dielectric.  A  cylindrical  base  that  has  a 
dielectric  constant  of  35  and  electric  conductivity  of  IS/m  was  added  to  the  dielectric 
hemisphere  to  dampen  the  reflections  from  the  boundary  which  is  opposite  to  the  antenna. 

The  modeling  results  show  that  the  pulsed  radiation  is  focused  at  the  second  focal  point,  which 
is  at  6cm  distance  from  the  hemisphere  surface  (Fig.4a).  Fig.  4b  shows  that  the  electric  field  at 
the  second  focus  (distance  =  6  cm)  is  50  V/m  for  a  ramp  input  of  1  V  for  a  zero  conductivity 
dielectric.  It  decreases  to  20  V/m  for  a  conductivity  of  0.3  S/m.  The  intensity  of  the  electric  field 
is  obviously  strongly  dependent  on  the  medium  conductivity:  higher  conductivity  causing 
stronger  attenuation  in  the  pulse  propagation  path  and  consequently  lower  field  intensity  at  the 
second  focus. 
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Fig.4.  (a)  A  1  V  ramp  voltage  with  a  risetime  of  100  ps  is  launched  and  focused  into  a  hemispherical 
target  with  diameter  of  12  cm.  The  focal  point  is  6  cm  from  the  sphere  surface;  b)  The  distribution  of 
lateral  electric  fields  in  the  axial  direction  (z-axis)  for  conductivities  from  0-0.5  S/m. 
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The  relatively  poor  matching  condition  on  the  air-dielectric  interface  causes  intense  reflection  of 
the  electromagnetic  waves  on  the  boundary,  reducing  the  intensity  of  the  transmitted  wave  even 
further.  Therefore,  it  is  important  to  place  a  matching  layer  with  a  permittivity  between  that  of  air 
and  the  dielectric  on  top  of  the  hemispherical  target  (or  any  other  target)  in  order  to  reduce  the 
reflections  at  this  interface. 

Increasing  the  focusing  resolution  using  a  dielectric  lens 

The  dielectric  properties  of  the  target  medium  and  its  surroundings  determine  to  a  large  extend 
how  effective  the  transfer  of  electromagnetic  energy  from  the  antenna  into  the  target  is.  If  the 
wave  incident  on  the  target  is  in  air,  but  the  target  medium  has  a  large  relative  permittivity  (say, 
about  81  for  water),  then  there  will  be  a  significant  reflection  of  the  pulse  (80%).  Adding  a 
dielectric  lens  composed  of  well  defined  dielectric  layers  to  the  antenna  allows  a  better  matching, 
and  consequently  a  higher  electric  field  in  the  target.  In  addition  it  allows  us  to  reduce  the  spot 
size,  thereby  increasing  the  field  on  the  target  even  further  (Fig. 5). 

The  lens,  which  was  designed  by  Carl  Baum,  is  hemispherical  [32-35]  and  the  measurement  in 
the  lens  dielectrics  was  discussed  in  [36-38],  Its  focal  point  is  in  the  geometric  center  of  the 
hemisphere  and  overlaps  the  focus  of  the  reflector  antenna.  The  lens  consists  of  five  layers  of 
dielectric  materials  with  dielectric  constants  varying  in  an  exponential  profile  from  free  space  to 
the  innermost  layer  srmax.  Due  to  the  increase  of  the  dielectric  constant,  the  focal  spot  size  can 
be  reduced  by  a  factor  of  8rmax("1/2)  at  the  innermost  layer.  The  choice  of  the  number  of  layers  and 
their  thickness  is  optimized  to  permit  maximal  transmission  through  these  layers.  The  innermost 
layer  is  generally  made  of  the  same  material  as  the  medium  to  be  probed.  In  this  lens,  the 
outermost  layer  has  a  dielectric  constant  of  1.3,  and  the  innermost  layer  has  a  dielectric 
constant  of  9.  The  spot  size  at  the  lens  focus  is  approximately  1/3  of  that  in  free  space. 

Fig. 6a  shows  the  increase  in  pulse  amplitude  right  at  the  surface  of  the  hemispherical  lens,  at 
the  point  z0.:  the  ratio  of  the  peak  electric  fields  at  this  point  compared  to  that  without  lens  in  air 
is  1.408[39].  Since  the  half-widths  of  the  impulse  responses  in  air  and  in  the  lens  are  not 
identical,  the  electric  enhancement  is  given  by  the  ratio  of  the  integral  of  the  curves,  which  can 
approximately  expressed  as  A=  (halfwidth)(peak  arnplitude)=FWHM*Emax.  According  to  this 
calculation  the  electric  enhancement  is  AimpU|Se.iens/Aimpuise-air®2.11.  Fig.  6b  compares  the 
numerical  simulation  and  experimental  results  of  the  beam  width  in  air  and  in  the  lens.  The 
beam  width  was  reduced  by  a  factor  of  two  from  4.45  cm  (no  lens)  to  2.04  cm  (with  lens). 


Fig. 5.  An  impulse  radiating  antenna  is  used  in 
combination  with  a  dielectric  lens  which  is  located  on 
the  surface. 


6 


Fig. 6.  a)  Comparison  of  impulse  waveforms,  at  the  second  focal  point,  in  air  and  at  the  lens  surface,  at  z0; 
b)  Comparison  of  simulation  and  experiment  results  of  beam  width  with  and  without  lens  (at  position  z0). 
The  pulse  risetime  is  -100  ps. 


Target  Detection  with  Focused  Subnanosecond  Radiations 

The  antenna  system  as  described  above  has  been  used  to  study  the  use  of  subnanosecond 
pulses  to  detect  abnormalities  in  a  homogeneous  background.  When  the  antenna’s  radiation  is 
focused  and  the  target  at  the  focus  is  scanned,  coherent  backscattered  signals,  signals  with  the 
largest  amplitude,  can  be  recorded  by  the  same  antenna.  The  incoherent  backscattered  signals 
from  the  out-of-focus  background  clutter,  are  shifted  in  time  and  can  be  separated  from  the 
target  signals.  The  technique  is  referred  to  as  confocal  imaging.  The  advantages  of  confocal 
imaging  are  increased  signal-to-noise  ratio  and  signal-to-clutter  ratio  and,  therefore,  improved 
target  detection  capability. 


Fig. 7  (a)  A  focusing  antenna  in  combination  with  a  dielectric  lens  illuminates  a  target  and  the 
backscattering  is  collected  with  the  same  antenna,  (b)  The  actual  setup. 
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In  the  confocal  imaging  system,  decreasing  the  focal  spot  size  is  important  to  increasing  the 
resolution.  One  way  to  reduce  the  spot  size  is  to  use  a  dielectric  lens  in  conjunction  with  the 
reflector  antenna,  as  shown  in  Fig.  7.  The  lens  is  hemispherical.  Its  focal  point  is  in  the 
geometric  center  of  the  hemisphere  and  overlaps  the  focus  of  the  reflector  antenna.  Therefore, 
the  antenna,  in  combination  with  the  lens,  constitutes  a  confocal  imaging  system. 

To  determine  the  imaging  resolution  of  the  antenna-lens  system,  we  have  imaged  two  targets  of 
equal  dimensions  2.5X1 .5X4. 5cm3  but  with  different  spacing.  The  two  objects  are  shifted  from 
one  edge  at  the  surface  of  the  lens  to  the  other.  This  is  equivalent  to  laterally  scanning  the  two 
objects.  In  the  experiment,  the  pulse  risetime  was  200  ps,  which  suggests  a  focal  spot  size  of  4 
cm,  derived  from  the  lOOps  data  shown  in  Fig. 6b. 

The  background  signals  are  subtracted  from  the  signals  obtained  at  each  target  location.  The 
value  of  the  difference  is  used  to  construct  a  2-D  image  by  assigning  pixel  values  according  to 
the  magnitude  of  the  signals  received.  Fig.  8  shows  the  image  constructed  with  two  objects  with 
spacing  of  3.5  and  5.5  cm  [40].  The  two  targets  are  irresolvable  in  the  lateral  direction  when  the 
spacing  is  less  than  3.5  cm.  However,  the  imaging  system  was  able  to  distinguish  targets  for  a 
spacing  of  5.5  cm.  The  resolution  was  therefore  estimated  to  be  approximately  4  cm. 

Our  present  confocal  system  prevents  a  scanning  of  targets  at  various  distances  from  the  lens. 
The  targets  need  to  be  placed  close  to  the  lens,  within  2  cm  of  the  focus,  in  order  to  provide 
return  signals  which  can  be  distinguished  from  those  obtained  without  lens.  Scanning  of  targets 
at  various  depths  may  be  possible  if  the  probed  object  is  placed  within  a  hollowed-out  lens. 
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Fig. 8.  Images  of  two  identical  metallic  objects  that  are  separated  by  (a)  3.5  cm  and  (b)  5.5  cm.  The  two 
objects  can  be  spatially  resolved  at  a  lateral  spacing  of  5.5  cm,  but  are  irresolvable  for  3.5  cm. 


8 


Modifying.  Target  Physical  Structure  with  Intense  Subnanosecond_  Pulses 
Focused  subnanosecond  pulsed  radiation  not  only  has  the  potential  for  detecting  an 
inhomogeneity  in  a  complex  background,  it  also  has  the  ability  to  modify  structures  of  properties 
of  targeted  objects,  such  as  biological  cells,  if  the  intensity  of  the  pulses  is  sufficiently  enhanced. 

It  is  known  that  for  pulses  with  durations  shorter  than  the  charging  time  of  the  outer  membrane 
(-100  ns  for  mammalian  cell),  the  probability  of  electric  field  acting  upon  intracellular  structure 
increases.  Shorter  pulses  in  the  subnanosecond  range  are  consequently  even  more  effective 
than  nanosecond  pulses  in  reaching  intracelluar  structures.  However,  this  proposition  should  not 
create  the  impression  that  cell  membrane  will  be  exempt  from  the  effect  of  subnanosecond 
pulses.  We  have  found  that,  for  the  first  time,  cell  membranes  increase  their  conductance,  in 
response  to  high  power  subnanosecond  electric  pulses  (20  kV/cm,  200  ps). 

Shown  in  Fig. 9  are  the  voltage-current  characteristics  of  a  cell’s  membrane  (NG108, 
Neuroblastoma  Cell  Lines)  after  it  is  exposed  to  subnanosecond  pulses  [41].  The  IV  curves 
were  obtained  by  means  of  the  patch  clamp  technology.  The  membrane  currents  were  recorded 
by  sweeping  the  voltages  from  -100  mV  to  40  mV.  At  10  and  30  seconds  after  exposure  of  2000, 
200  ps,  20  kV/cm  pulses  delivered  at  10  Hz,  the  cell  membrane  becomes  more  permeable:  it 
shows  a  higher  leakage  current  compared  with  the  cell  that  is  not  pulsed.  The  leakage  current 
only  flows  in  the  inward  direction  from  the  extracellular  to  intracellular  medium.  The  pathway  of 
the  current  flow  is  most  likely  to  be  the  small  pores  formed  in  the  cell  membranes,  as  is  proved 
through  other  membrane  integrity  markers,  TI+  (atomic  diameter:  0.34  nm)  and  Propidium+ 
(1.38x1.15  x0.54  nm).  These  pores  were  estimated  on  the  order  of  -1.4  nm  in  diameter. 
Unexpectedly,  the  current  rectifying  characteristic  of  these  pores  is  similar  to  synthetic 
nanopores  in  polymer  foils,  reportedly  to  be  in  asymmetrical  conical  shape  [42], 

Whereas  causing  cell  membrane  leakage  requires  2000  pulses,  we  can  also  induce  cell  death 
by  increasing  the  number  of  pulses  with  a  high  repetition  rate  (10  kHz  in  burst  mode).  We  have 
observed  that  subnanosecond  pulses  caused  cell  death  with  a  significant  temperature  rise  due 
to  a  high-repetition  rate  and  high  pulse  numbers.  Cell  viability  assessed  by  WST  assay  18  h 
after  the  exposure  of  2.4  million  pulses  was  found  to  be  34%+  2%  with  temperature  rising  to 
42°C^45°C  during  exposure. 

The  cell  death  is  due  to  both  electric-field  and  thermal  effects.  Delivering  pulses  of  higher 
electric  field  intensity  at  a  lower  repetition  rate  may  allow  us  to  separate  these  two  effects,  which 
justifies  further  study.  On  the  other  hand,  subnanosecond  pulses  offer  a  potential  treatment 
modality  that  can  modify  cell  membrane  properties  or  alter  the  physiological  conditions  of  cells 
to  induce  a  desired  response  instead  of  causing  cell  death.  In  this  case,  much  lower  numbers  of 
pulses  are  needed  to  obtain  this  effect.  Such  applications  can  include  enhancement  of  wound 
healing,  or  hemostasis  (blood  clotting),  by  stimulating  calcium  influx  into  platelets. 
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Leak-current 


Fig. 9.  Current-voltage  characteristics  of  a  cell’s 
membrane  after  being  exposure  to  200  ps,  20  kV/cm 
pulses.  According  to  the  convention  in  the  patch 
clamp,  the  positive  current  flows  outward  from  inside 
to  outside  cell  (in  the  first  quadrant).  Conversely,  the 
negative  current  flows  inward  (in  the  third  quadrant) 
from  outside  to  inside  cell.  After  exposure  of  2000 
pulses,  the  cell  shows  leakage  current  increase  in 
negative  direction  but  no  change  in  positive  direction, 
which  indicates  inward  “diodes”  are  produced  in  the 
cell  membrane. 


Future  Work 

Improving  the  spatial  resolution  for  imaging 
targets  and/or  inducing  physical  changes  in 
targets  is  possible  but  requires  modifications  of  our  system:  The  target  and  its  surrounding 
dielectrics  would  need  to  be  included  into  the  imaging  system  as  part  of  the  lens.  For  inducing 
the  change  of  dielectrics,  a  very  high  power  wave  launcher  at  the  first  focus  needs  to  be 
constructed  [43-44],  The  material  of  the  focusing  lens  has  a  tapered  profile  of  dielectric 
constants  in  order  to  match  the  surrounding  dielectrics.  This  lens  will  still  be  fed  with  the 
reflector  antenna,  but  it  would  need  to  be  hollowed  out  near  the  focus,  by  creating  a 
hemispherical  or  spheroidal  cavity  (Fig.  10).  The  dielectric  medium  surrounding  the  target  would 
then  be  embedded  into  this  cavity  [45]. 

Using  this  technique,  allows  us  to  increase  the  solid  angle  through  which  the  electromagnetic 
wave  can  approach  the  target  through  the  lens.  If  the  dielectric  medium  surrounding  the  target  is 
pliant,  it  could  then  be  sucked  into  this  cavity  by  several  small-diameter  tubes  drilled  through  the 
lens  media  to  the  lens/dielectrics  interface.  Transverse  stretching  will  reduce  the  dispersion  of 
the  pulse  as  it  approaches  the  target.  This  technique  may  be  limited  to  not-too-deeply-buried 
targets.  As  one  goes  to  deeply  buried  targets,  one  will  need  to  increase  the  radius  of  the  lens 
aperture  to  sufficiently  large  values  compared  to  the  depth.  This  will,  in  turn,  increase  the  overall 
lens  radius. 
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Our  ultimate  goal  is  to  develop  a  useful  delivery  method  of  subnanosecond  pulses  that  is  able  to 
deposit  a  substantial  amount  of  energy  or  monitor  a  specific  area  at  a  specific  depth  inside  a 
dielectric  subject  placed  at  the  prolate  spheroid  reflector’s  focal  point. 


Fig.  10.  Embedding  dielectric,  hemispherical  target  such  that  it 
becomes  part  of  the  lens  system  allows  us  to  increase  the  solid 
angle  through  which  the  electromagnetic  wave  can  approach 
the  target  through  the  lens. 
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Abstract — A  prolate-spheroidal  reflector  antenna  focusing  sub¬ 
nanosecond  pulsed  electromagnetic  radiation  in  the  near  field  is 
studied.  This  antenna  reflector  is  fed  by  a  pair  of  balanced,  con¬ 
ical-plate  transmission  lines  and  the  radiated  pulse  consists  of  a 
prepulse,  an  impulse,  and  a  postpulse.  Whereas  the  prepulse  am¬ 
plitude  decreases  inversely  with  distance,  the  impulse  is  maximum 
near  the  geometric  focus,  and  its  power  density  distribution  has  a 
full- width  at  half-maximum  (FWHM)  width  of  32  cm  in  the  axial 
direction  and  10  cm  in  the  lateral  direction. 

Index  Terms — Focusing,  impulse  radiating  antennas  (IRAs), 
near  field,  prolate  spheroidal  reflector,  subnanosecond  pulse. 


I.  Introduction 

RESEARCH  on  the  near-field  radiation  of  subnanosecond 
electric  pulses  is  motivated  by  the  possible  use  of  sub¬ 
nanosecond  pulsed  antennas  for  medical  therapies  and  dielectric 
imaging.  Biological  cells  exposed  to  the  subnanosecond  electric 
pulses  have  shown  increased  lethality  and  increased  calcium  re¬ 
lease,  which  may  lead  to  medical  therapies,  such  as  cancer  treat¬ 
ment  or  wound  healing  [1] .  At  lower  intensity,  antennas  fed  with 
subnanosecond  pulses  can  be  used  to  explore  the  electrical  prop¬ 
erties  of  biological  tissues  in  the  near-field  region,  allowing  us 
to  obtain  the  dielectric  profile  of  the  tissue  under  observation 
and  to  detect  abnormalities  in  a  uniform  background  due  to  the 
dielectric  contrast  [2],  [3].  Both  imaging  and  medical  applica¬ 
tion  require  focusing  the  radiation  to  a  minimum  spot  size.  In 
imaging,  a  small  spot  size  allows  for  a  3-D  scan  with  a  good 
spatial  resolution.  For  therapy,  a  small  spot  size  means  that  only 
the  area  at  the  focal  point  will  be  treated  as  opposed  to  the  sur¬ 
rounding  area  with  lower  electric  field  (or  power  density). 

The  antennas  for  radiating  and  focusing  ultrawideband, 
high-power  signals  with  low  dispersion  and  high  directivity 
can  be  reflector-type  antennas,  lenses  or  antenna  arrays.  For 
reflector-type  antennas,  an  ellipsoidal  dish  is  usually  used. 
The  dish  has  two  focal  points  with  a  source  placed  at  the  first 
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focal  point  and  the  target  placed  at  the  second  focal  point.  In 
research  on  microwave-induced  hyperthermia,  an  ellipsoidal 
reflector  for  focusing  a  narrowband  signal  was  used  [4].  For 
high-power  ultrawideband  radiation  in  the  time  domain,  studies 
were  mostly  focused  on  long-range,  far-field  impulse  radiating 
antennas  (IRAs),  in  which  a  paraboloidal  dish  with  a  second 
focus  at  infinity  is  employed  [5].  For  the  second  focal  point 
close  to  the  near-field  region  with  focal  distance  less  than 
2D2 /A  (D:  the  aperture  diameter;  A:  the  wavelength),  the 
spot  size  is  reduced,  and  the  spot  width  W  is  approximately 
proportional  to  the  focal  distance  (ro)  :  W  =  roA/D,  which 
indicates  that  a  large  aperture  for  a  short  focal  distance  permits 
a  small  focus  [6],  [7].  Recently,  a  series  of  studies,  using  the 
prolate  spheroidal  reflectors  in  conjunction  with  the  typical 
TEM  feed-arms  used  in  long-range  IRA,  have  shown  focusing 
of  a  radiated  impulse  in  the  near  field  [8],  [9].  As  reported  in  [8] 
and  [9],  the  antenna  had  a  1-m  aperture  diameter  with  a  quarter 
of  the  prolate  spheroid  as  a  reflector.  The  reflector  was  installed 
on  a  ground  plane  and  an  unbalanced  feed  structure  was  used. 

In  this  letter,  we  report  a  smaller  version  of  a  prolate  sphe¬ 
roidal  IRA  (PSIRA)  with  an  aperture  diameter  (0.5  m)  designed 
for  near-field  imaging  or  for  therapeutic  applications.  Whereas 
the  previously  reported  PSIRA  was  primarily  studied  for  the 
proof  of  concept,  this  PSIRA  was  designed  a  step  closer  to  the 
actual  applications,  which  require  the  antenna  to  be  compact 
and  portable.  To  our  knowledge,  this  is  the  first  IRA  with  a  bal¬ 
anced  feed  configuration,  which  focuses  the  impulse  in  the  near 
field.  This  antenna  uses  half  of  the  prolate  spheroid  as  the  re¬ 
flector.  The  antenna  structure  was  scaled  down  by  a  factor  of 
0.5  from  the  suggested  PSIRA  dimensions  [8].  It  was  manufac¬ 
tured  by  Farr  Research,  Albuquerque,  NM.  We  have  measured 
the  electric  field  distribution  of  this  antenna  with  an  approxi¬ 
mate  Gaussian  input  of  half- width  200  ps. 

II.  Antenna  Structure 

The  antenna  consists  of  three  parts:  a  prolate  spheroidal  re¬ 
flector,  a  conical  feed  structure,  and  a  balun  (Fig.  1).  The  semi¬ 
major  axis  of  the  reflector,  a  =  29.8  cm,  and  the  semiminor  axis 
b  =  25  cm.  The  focal  distance  zq  is  16.3  cm.  The  conical  feed 
has  two  pairs  of  parallel,  balanced  conical  plates  as  a  transmis¬ 
sion  line,  each  arranged  at  an  angle  of  60°  with  respect  to  the 
y-axis  resulting  in  an  impedance  of  200  £1.  A  balun,  made 
of  two  100-f)  cables,  allows  a  1:4  ratio  impedance  transforma¬ 
tion  from  the  pulse  input  (50  ft  from  the  pulse  source)  to  the 
conical  transmission  line  (200  f)).  The  balun  has  the  advantage 
of  doubling  the  voltage  at  the  input  terminal  of  the  transmission 
line.  The  field  at  the  focal  point,  with  our  feed-arms  at  60°,  is 
amplified  by  1.6  compared  to  one  pair  of  conical  feed-arms  [9]. 
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Fig.  2.  The  schematics  of  the  analytical  pulse  shape  at  the  focus  when  the  an¬ 
tenna  is  fed  with  a  step-function  pulse. 


on  the  surface  will  arrive  at  the  second  focal  point  in  the  same 
phase. 

Shown  in  Fig.  2  is  the  analytical  focal  waveform  expected 
at  the  focus  for  a  fast-rising  step  function  feeding  the  antenna. 
The  impulse  part  primarily  contains  E s,  subtracted  by  Ep.  Es 
is  proportional  to  the  time  derivative  of  the  feed  pulse:  A  faster 
rise-time  pulse  produces  a  higher  amplitude.  The  prepulse,  Ep? 
is  proportional  to  the  amplitude  of  the  feed  pulse.  The  step  im¬ 
pulse,  Es,  from  the  conical  transmission  in  the  opposite  direc¬ 
tion  to  Ep,  is  reflected  on  the  reflector  and  also  contributes  to 
the  impulse.  Es  and  Es  are  in  the  same  direction.  Ep  is  valid 
until  the  wave  travels  to  the  truncation  of  the  reflector  along  the 
conical  transmission  line.  After  that,  Ep  becomes  a  new  term, 
EP2.  The  analytical  analysis  of  these  terms  can  be  seen  in  [9]. 
In  the  impulse,  we  are  interested  in  the  dominant  part,  E^.  For  a 
given  step-input,  the  x-axis-directed  (lateral)  components  of  the 
prepulse  and  the  impulse  at  the  geometric  focus  (zo  =  16  cm) 
were  calculated  by  integrating  the  fields  of  the  reflected  TEM 
wave  over  the  truncated  aperture  plane  [8] 

Ep  =  1.6[V(t)/(27rfg)][(a  + z0)/(bz0)]  (1) 

Es  *  (1.6/c)  [d V (t ) /dt]  [(b/a) / (irfg )] .  (2) 


Fig.  1 .  Antenna  structure,  (a)  Side  view.  A  differential  signal  is  fed  to  the  trans¬ 
mission  line  at  the  first  focal  point  (— z0).  (b)  Front  view.  ©0  is  the  angle  be¬ 
tween  the  feed-arm  and  y-axis.  (c)  The  actual  antenna. 


The  waves  scattered  on  the  reflector  surface  become  an  exact 
inhomogeneous  spherical  TEM  wave  propagating  toward  the 
second  focus.  At  the  end  of  the  conical  feed  plate,  resistors  were 
placed  to  absorb  the  low-frequency  component  of  the  signals. 

III.  Electric  Field  at  the  Focal  Point 

The  electric  field  at  the  second  focal  point  comprises  three 
components:  the  prepulse,  the  impulse,  and  the  postpulse.  The 
prepulse  (Ep)  is  radiated  directly  from  the  conical  transmission 
line,  and  the  impulse  is  the  electromagnetic  pulse  reflected  from 
the  prolate- spheroidal  reflector.  Only  the  impulse  will  be  fo¬ 
cused  at  the  second  focal  point  coherently  as  the  reflected  wave 


V(t)  is  the  voltage  of  the  feeding  pulse,  c  is  the  speed  of  light, 
and  fg  is  the  characteristic  impedance  factor  of  the  conical  feed- 
arms.  It  is  noted  that  electric  fields  are  multiplied  by  a  factor  of 
1.6  for  the  antenna  with  four-feed  arms  at  60°.  For  the  impulse 
part,  the  duration  of  the  impulse  is  t  and  the  field  amplitude  of 
the  impulse  is  proportional  to  Vot^  1  (V/m)  for  a  step  function 
with  an  amplitude  of  Vo  and  a  rise  time  of  t$. 

From  (2),  it  is  clear  that  in  order  to  produce  a  large  field  at  the 
focal  point,  a  large  voltage-time  derivative  and  a  low  launching 
impedance  of  the  conical  transmission  line  (low  fg)  are  required. 
The  pulse  rise  time  should  be  short  in  order  to  increase  the  am¬ 
plitude  of  the  impulse. 

The  impulse  and  the  prepulse  have  opposite  signs,  with  the 
impulse  being  much  larger  than  the  prepulse  (Fig.  2).  The  ratio 
of  the  impulse  to  the  prepulse  is  zo  T/a,  in  which  T  satisfies 
T  =  2 [a  —  z0]/ct6.  This  ratio  suggests  that  in  order  to  obtain 
a  large  impulse,  zo/a  needs  to  be  large.  A  “fat”  and  shallow 
antenna  should  produce  a  larger  impulse  amplitude  at  the  focus. 
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Fig.  3.  (a)  The  measured  electric  pulse  from  the  pulse  generator  that  was  fed 
into  the  antenna,  (b)  The  measured  electric  field  at  the  second  focal  point  (solid 
line).  The  dash  line  is  dV/dt,  and  V  is  the  feed  pulse  shown  in  (a).  The  peak  of 
dV/dt  was  normalized  to  the  impulse  amplitude. 


IV.  Electric  Field  Distribution 

To  measure  the  distribution  of  the  electric  field  at  and  around 
the  focal  point  of  the  PSIRA,  we  used  a  TEM  horn  antenna 
with  two  tapered,  curved  plates,  which  has  the  same  polariza¬ 
tion  as  the  PSIRA  and  measures  the  x-directed  electric  field. 
The  antenna  has  an  opening  distance  of  4  cm.  It  was  facing  to¬ 
ward  the  prolate  spheroidal  antenna  aperture.  The  antenna  was 
moved  along  the  z-axis  to  measure  the  longitudinal  electric  field 
distribution.  Fig.  3(a)  shows  the  input  pulse,  generated  by  a 
subnanosecond  pulse  generator  (FPG5-10PM,  manufactured  by 
FID).  As  the  sensor  measures  the  time  rate  change  of  electric 
displacement  field  [10],  the  measured  signal  was  integrated  in 
time  to  obtain  the  electric  field  [Fig.  3(b)]. 

In  Fig.  3(b),  the  prepulse  is  wider  than  the  impulse  since  the 
prepulse  is  just  the  replica  of  the  input  pulse  and  the  impulse 
is  proportional  to  the  time  derivative  of  the  feed  pulse  shown 
in  Fig.  3(a).  The  feed  pulse,  which  was  directly  measured  from 
the  pulse  generator,  was  differentiated  to  obtain  dV/dt  [dashed 
line  in  Fig.  3(b)].  Apparently,  the  impulse  part  of  the  electric 
field  is  overlapped  with  dV/dt,  which  confirms  the  prediction 
of  (2).  The  difference  of  frequency  components  of  prepulse  and 
impulse  allows  us  to  use  a  high-pass  filter  to  suppress  the  pre¬ 
pulse  and  enhance  the  impulse  in  future  applications,  such  as 
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Fig.  4.  (a)  The  measured  electric  field  (x-directed)  distribution  in  the  axial  di¬ 
rection  (z-axis).  The  prepulse,  the  impulse,  and  the  square  of  impulse  are  shown. 
The  curve  with  down-triangle  symbols  is  the  modeled  impulse  distribution  in 
z-axis.  (b)  The  measured  electric  field  (x-directed)  distribution  in  the  lateral  di¬ 
rection  (y-axis).  The  prepulse,  the  impulse,  and  the  square  of  impulse  are  shown. 
The  curve  with  down-triangle  symbols  is  the  modeled  impulse  distribution  in 
y-axis. 


in  a  near-field  radar.  The  separation  in  time  between  these  two 
pulses,  defined  as  cleartime,  is  approximately  1  ns.  The  post¬ 
pulse  after  t  =  4.5  ns  may  contain  both  diffractions  from  var¬ 
ious  parts  of  the  antenna  (reflector  rim,  edge  of  the  conical  plate, 
etc.)  and  the  resonance  signal  of  the  sensor  itself. 

The  amplitudes  of  both  prepulse  and  impulse  are  plotted  in  y- 
and  z-axes  in  Fig.  4.  The  prepulse  decreases  in  the  z-direction, 
with  a  1/z  dependence.  For  an  observation  point  at  z  =  16  cm, 
the  geometric  focus,  the  prepulse  stays  approximately  constant 
along  the  y-axis  (y  =  —12  to  12  cm),  indicating  again  that 
it  scales  as  1/z.  The  power  density  of  the  impulse,  which  is 
proportional  to  the  square  of  the  electric  field,  is  focused  in  both 
lateral  and  axial  directions.  The  distribution  of  power  density 
has  a  full- width  at  half-maximum  (FWHM)  width  of  32  cm  in 
the  axial  direction  and  a  FWHM  width  of  10  cm  in  the  lateral 
direction  (at  z  =  16  cm). 

The  electric  field  distribution  of  this  antenna  was  also  mod¬ 
eled  by  CST  Microwave  Studio,  a  3-D,  finite  integral  time  do¬ 
main  (FITD)  electromagnetic  simulation  software.  The  input 
excitation  is  a  Gaussian  waveform  with  a  half-width  of  200  ps  to 
closely  model  the  experimental  input  waveform.  The  modeled 
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distributions  of  the  impulse  in  the  axial  direction  and  the  lateral 
direction  are  plotted  in  Fig.  4(a)  and  (b),  respectively.  The  mod¬ 
eled  results  agree  very  well  with  the  experimental  results.  This  is 
true  especially  near  the  focal  point  in  both  the  axial  and  the  lat¬ 
eral  electric  field  distributions.  At  an  observing  distance  closer 
to  the  reflector  edge,  we  observe  a  slight  discrepancy  between 
the  modeled  field  and  measured  field,  which  could  be  due  to  the 
fact  that  the  sensor  may  pick  up  other  electric  field  components 
besides  the  x-directed  field  and  the  model  only  gives  the  x-di- 
rected  field. 

We  also  note  that  the  maximum  of  the  impulse  electric  field 
on  the  axis  is  shifted  slightly  from  the  geometric  focus  (z  = 
16  cm)  toward  the  reflector,  which  is  observed  in  both  the  mea¬ 
sured  result  and  the  modeled  result.  This  is  because  the  im¬ 
pulse  decreases  inversely  with  the  distance  while  it  is  focused 
in  space.  At  the  focal  point,  even  though  we  have  a  coherent 
combination  of  waves,  the  impulse  electric  field  is  still  smaller 
than  the  nearby  locations  toward  the  reflector  due  to  a  large  im¬ 
pulse  width,  ct<5  (in  spatial  units).  In  order  that  the  maximum 
impulse  amplitude  occurs  at  the  geometric  focus,  the  impulse 
width  needs  to  be  small  compared  to  both  2zo  and  2b  (similar 
discussion  in  the  frequency  domain  can  be  seen  in  [7]).  A  pulse 
with  faster  rise  time  should  allow  the  shift  of  the  focal  spot  to¬ 
wards  the  geometric  focus. 

V.  Conclusion  and  Future  Work 

The  PSIRA  with  a  reflector  diameter  of  0.5  m  focuses  the  sub¬ 
nanosecond  electric  pulse  radiation  near  the  second  focal  point 
(z  =  16  cm).  The  focus  spot  has  an  ellipsoidal  shape,  and  the 
distribution  of  power  density  has  a  small  diameter  in  the  lat¬ 
eral  direction  and  a  wide  dimension  in  the  axial  direction  (10 
and  32  cm  in  FWHM).  The  measured  electric  field  distribution 
agrees  with  the  modeled  results  to  a  large  extent.  Near  the  focal 
point,  the  measured  electric  field  is  exactly  as  the  model  pre¬ 
dicts.  When  the  focusing  antenna  is  used  for  the  delivery  of 


short  pulses  into  targets  in  the  near  field,  a  focusing  dielectric 
lens  near  the  second  focal  point  may  be  used  to  minimize  the 
spot  size  and  increase  the  impulse  amplitude  [11]. 
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Target  Detection  With  Impulse  Radiating  Antenna 

Chandra  Bajracharya,  Shu  Xiao,  Carl  E.  Baum,  and  Karl  H.  Schoenbach 


Abstract — A  prolate  spheroidal  reflector  antenna  focuses 
subnanosecond  pulsed  radiation  in  the  near  field  with  a  small 
beamwidth,  allowing  for  applications  such  as  target  detection. 
Confocal  imaging  in  free  space  has  been  demonstrated  with  such 
an  antenna  for  both  transmitting  and  receiving.  A  polarized  beam 
allows  for  the  detection  of  linear  target  orientation.  We  show  that 
decreasing  the  focal  spot  size  by  adding  a  dielectric  lens  enables 
us  to  resolve  a  target  with  a  typical  dimension  of  4  cm.  The  spatial 
resolution  can  be  further  enhanced  by  increasing  the  dielectric 
constant  of  the  lens’  innermost  layer. 


Antenna 


Index  Terms — Confocal  imaging,  dielectric  lens,  near  field, 
subnanosecond  pulses. 


Oscilloscope 


I.  INTRODUCTION 


MICROWAVE  imaging,  either  with  narrowband  or  wide¬ 
band  pulses,  can  detect  the  presence  of  abnormalities 
in  a  dielectric  medium  or  retrieve  the  dielectric  properties  of  a 
target  through  an  inverse  scattering  technique.  Wideband  pulses 
in  imaging  systems  have  been  used  in  subsurface  probing,  med¬ 
ical  imaging,  and  seeing  through  walls  [1]— [3]  because  of  the 
high  resolution  (centimeters  in  free  space)  and  the  ability  to  pen¬ 
etrate  obstacles  over  a  broad  spectrum  of  electromagnetic  radi¬ 
ation.  Scanning  a  region  of  interest  can  be  performed  by  one 
antenna  using  synthetic  aperture  methods  [4]  or  by  an  antenna 
array  through  digital  beamforming  [5].  When  the  antenna’s  ra¬ 
diation  is  focused  and  the  target  at  the  focus  is  scanned,  coherent 
backscattered  signals,  the  signals  with  the  largest  amplitude, 
will  be  recorded  by  the  same  antenna.  The  incoherent  backscat¬ 
tered  signals  from  the  out-of-focus  background  clutter,  on  the 
other  hand,  are  shifted  in  time  and  can  be  separated  from  the 
target  signals.  The  technique  is  referred  to  as  confocal  imaging. 
The  advantages  of  confocal  imaging  are  increased  signal-to- 
noise  ratio  and  signal-to-clutter  ratio  and,  therefore,  improved 
target  detection  capability.  The  concept  of  confocal  imaging 
using  a  fixed- focus  antenna  was  discussed  in  [2],  but  the  details 
of  the  design  were  not  discussed,  primarily  because  the  interest 
was  focused  on  the  antenna  array,  which  has  the  ability  to  vary 
the  focus  in  using  the  digital  beamforming. 
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Fig.  1.  Confocal  imaging  system  with  a  near-field,  impulse  focusing  antenna. 
The  pulse  generator  delivers  pulses  with  a  rise-time  of  approximately  200  ps  to 
the  antenna,  which  focuses  the  radiation  in  a  focal  plane  16.5  cm  away  from  the 
aperture  plane.  The  prolate  spheroidal  reflector  has  a  major  axis  of  30  cm  and 
a  minor  axis  of  25  cm.  The  aperture  plane  with  a  diameter  of  50  cm  is  at  the 
mid-plane  of  the  reflector.  The  signals  scattered  by  the  target  in  free  space  are 
collected  by  the  same  antenna  and  sent  to  the  oscilloscope  for  analysis. 

In  this  letter,  we  report  the  use  of  a  fixed-focus  impulse  ra¬ 
diating  antenna  [6],  [7]  for  the  confocal  target  detection.  Target 
detection  using  this  system  is  limited  to  targets  in  free  space. 
Another  possible  use  would  be  to  detect  a  target  embedded  in  a 
weakly  scattering  object  so  the  focal  point  can  be  predetermined 
as  in  free  space.  We  also  use  a  lens  in  conjunction  with  the  an¬ 
tenna  to  decrease  the  focal  spot  size  and  to  increase  the  resolving 
ability.  This  lens-antenna  system  applies  in  the  case  when  a 
target  is  in  the  subsurface  region  of  a  medium  with  known  di¬ 
electric  properties.  The  dielectric  lens,  in  conjunction  with  the 
antenna,  can  be  used  to  match  the  impedance  from  free  space 
to  the  medium  and  to  achieve  the  focus  inside  the  medium.  One 
scenario  would  be  detecting  or  monitoring  cancers  in  the  skin 
with  a  lateral  scanning  at  a  fixed  depth. 

II.  Confocal  Imaging  System 

The  confocal  imaging  system  consists  of  a  pulse  generator 
(FPG-5P,  manufactured  by  FID  GmbH,  Germany),  an  oscillo¬ 
scope  (Tektronix  TDS7404),  and  a  reflector  antenna  (Fig.  1). 
The  reflector  antenna  is  based  on  two  foci  of  an  ellipse  and 
launches  the  subnanosecond  pulse  from  the  first  focal  point  to 
the  second,  which  is  16.5  cm  away  from  the  aperture  plane  [7]. 
In  the  confocal  configuration,  the  antenna  works  as  both  trans¬ 
mitter  and  receiver.  Usually,  the  operation  of  a  single  antenna  as 
transmitter  and  receiver  requires  an  ultrafast  transmit-receive 
switch  to  provide  adequate  isolation  between  transmit  and  re¬ 
ceive  ports.  Here,  we  instead  use  a  T-connector,  which  allows 
the  pulses  to  be  fed  to  the  antenna  and  also  allows  the  returned 
signal  to  be  sent  to  an  oscilloscope.  A  high-power  attenuator 
(NMFP-26B,  Bath  Electronics)  is  used  to  reduce  the  intensity 
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Fig.  2.  Recorded  signal  on  the  real-time  oscilloscope.  The  signals  consist  of 
four  parts:  1)  the  input  pulse  from  the  pulse  generator  (14%-92%  rise-time, 
200  ps,  FWHM:  ~600  ps);  2)  the  reflection  from  the  pulse  generator  because 
of  the  T-connector;  3)  the  reflected  signal  from  the  antenna  balun;  and  4)  the 
backscattered  signals  from  the  target  at  the  focus. 


of  the  signal,  followed  by  two  additional  wideband  attenuators 
(PE7045-20,  Pasternack)  on  the  scope. 

To  demonstrate  the  antenna’s  operation  as  a  confocal  trans¬ 
mitter  and  receiver,  we  place  an  aluminum  rod  (diameter: 
2.54  cm;  length:  15.5  cm)  at  the  focus  and  record  the  signals 
from  the  target.  Fig.  2  shows  the  signals  recorded  by  the 
oscilloscope  with  and  without  the  presence  of  the  target.  The 
recorded  signal  consists  of  four  main  parts:  the  input  pulse 
directly  from  the  pulse  generator,  the  reflections  from  the  pulse 
generator,  the  reflected  signal  from  the  antenna  balun,  and  the 
backscattered  signals  from  the  target.  The  first  three  signals 
should  be  invariant  and  independent  of  the  target  scattering. 
This  is  confirmed  by  overlapping  the  first  three  signals  for  both 
cases  with  and  without  the  target.  The  fourth  part  of  the  signal 
varies  as  a  result  of  the  introduction  of  the  target  and  contains 
backscattered  signals  from  the  target.  A  more  detailed  analysis 
of  the  signals  is  not  necessary  since  we  are  only  interested  in 
examining  the  changes  of  the  signals  in  the  presence  of  the 
target  versus  without  the  target. 

The  appearance  of  the  signal  changes  indicates  the  presence 
of  a  target  at  the  focal  point.  Once  the  time  window  of  this 
change  is  identified,  the  oscilloscope  is  set  in  a  delayed  mode 
so  that  this  particular  time  window  is  enlarged  by  setting  a 
lower  voltage  scale  and  shorter  timescale.  The  figure  inserted  in 
Fig.  2  shows  a  significant  difference  between  the  signals  with 
and  without  the  presence  of  the  target.  We  note  that  the  signal 
near  43-44  ns  for  the  case  without  a  target  is  most  likely  due 
to  the  reverberation  between  the  TEM  feed  and  the  reflector. 
As  a  result,  the  received  signal  by  the  same  antenna  contains 
ripples  in  the  tail,  which  coincide  with  the  time  window  of  the 
signal  returned  from  the  target.  This  coincidence  may  pose  a 
limitation  to  some  sensing  applications  where  the  change  of 
signal  obtained  by  subtracting  the  background  signals  is  not 
adequate.  One  way  to  reduce  the  ripples  is  to  employ  an  offset 
feed  arm  to  replace  the  centered  feed  arm,  as  discussed  in  a 
parabolic  impulse  antenna  [8]. 

One  attractive  feature  of  the  antenna  is  that  the  wave  at  the 
second  focus  is  linearly  polarized  and  can  be  adjusted  either  in 
the  vertical  or  horizontal  direction  by  rotating  the  antenna.  For 


Time  (ns) 

(a) 


Time  (ns) 

(b) 

Fig.  3.  Received  signal  for  different  orientations  of  a  target  when  the  antenna 
is  vertically  polarized,  (a)  The  rod  is  placed  at  the  focus  vertically  (copolarized), 
(b)  The  rod  is  placed  at  the  focus  horizontally  (cross-polarized). 


linear  targets,  e.g.,  a  metallic  rod  (diameter:  2.54  cm;  length: 
15.5  cm),  the  backscattered  signal  is  maximum  when  the  target 
is  copolarized  [Fig.  3(a)].  For  linear  targets  that  are  cross-po¬ 
larized,  the  antenna  rejects  the  backscattered  signals  [Fig.  3(b)], 
so  the  target  is  undetectable  to  the  imaging  system.  This  prop¬ 
erty  allows  us  to  selectively  scan  the  targets  by  examining  the 
intensity  of  the  target  return.  The  orientation  of  the  antenna  can 
be  varied  until  a  maximum  backscattered  signal  is  observed  and 
the  antenna  polarization  is  aligned  with  the  target  orientation. 

III.  Confocal  Imaging  System  With  a  Dielectric  Fens 

In  the  confocal  imaging  system,  decreasing  the  focal  spot  size 
is  instrumental  to  increasing  the  resolution.  One  way  to  reduce 
the  spot  size  is  to  use  a  dielectric  lens  in  conjunction  with  the  re¬ 
flector  antenna,  as  shown  in  Fig.  4.  The  lens  is  hemispherical.  Its 
focal  point  is  in  the  geometric  center  of  the  hemisphere  and  over¬ 
laps  the  focus  of  the  reflector  antenna.  Therefore,  the  antenna,  in 
conjunction  with  the  lens,  is  still  a  confocal  imaging  system.  The 
lens  consists  of  five  layers  of  different  dielectric  materials  with 
dielectric  constants  varying  in  an  exponential  profile  from  free 
space  to  the  innermost  layer,  £rmax  [9] .  Due  to  the  increase  of  the 
dielectric  constant,  the  focal  spot  size  can  be  reduced  by  a  factor 
of  at  the  innermost  layer.  The  choice  of  the  number  of 

layers  and  their  thickness  is  optimized  to  permit  maximal  trans¬ 
mission  through  these  layers.  The  innermost  layer  is  generally 
made  of  the  same  material  as  the  medium  to  be  probed.  In  this 
lens,  the  outermost  layer  has  a  dielectric  constant  of  1.3,  and 
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Fig.  4.  Reflector  antenna  used  in  conjunction  with  a  dielectric  lens.  The  focal 
point  of  the  dielectric  lens  at  the  center  of  the  hemisphere  is  at  the  same  position 
as  that  of  the  reflector  antenna.  The  diameter  of  the  lens  is  30  cm.  The  target  is 
shifted  along  the  lateral  axis  and  the  longitudinal  axis,  which  is  equivalent  to 
moving  the  antenna  and  lens  system  to  scan  the  space  along  the  front  of  the 
lens. 


Fig.  5.  Comparison  of  backscattered  signal  differences  with  and  without  the 
lens.  The  difference  is  obtained  from  the  backscattered  signal  after  subtraction 
of  the  background  signal.  The  target  is  a  metallic  object  (2 .5  x  1.5  x  4.5  cm3  ) , 
placed  at  the  focus.  Dashed,  horizontal  lines  indicate  the  peak  value  of  the 
difference. 


the  innermost  layer  has  a  dielectric  constant  of  9.  The  five-layer 
hemispherical  focusing  lens  was  made  to  the  lens  specifica¬ 
tions  [9]  by  TPL,  Inc.,  Albuquerque,  NM.  The  nanocomposite 
dielectric  materials  developed  by  TPL,  Inc.,  were  used  for  the 
fabrication  of  the  multilayer  lens  since  desired  dielectric  con¬ 
stants  can  be  achieved.  The  spot  size  at  the  lens  focus  is  approx¬ 
imately  1/3  of  that  in  free  space.  The  lateral  focal  spot  size  of 
the  reflector  antenna  was  measured  as  11  cm  [7],  so  at  the  lens 
surface  it  is  3-4  cm. 

The  antenna  and  lens  imaging  system  were  designed  to  work 
in  the  case  when  a  target  is  embedded  in  a  medium  that  has  the 
same  dielectric  constant  as  the  innermost  layer  [9] .  Therefore, 
this  imaging  system  requires  knowledge  of  the  medium  dielec¬ 
tric  property  when  designing  such  a  lens.  Here,  we  consider  a 
simpler  case:  a  metallic  target  (2.5  x  1.5  x  4.5  cm3)  placed  in 
free  space  and  in  front  of  the  lens.  The  experimental  setup  is 
shown  in  Fig.  4.  We  compare  the  difference  of  the  backscat¬ 
tered  signals  with  and  without  the  lens.  As  shown  in  Fig.  5, 
the  lens  significantly  increases  the  difference  of  the  scattered 
signal.  The  maximum  difference  (peak  to  peak)  for  the  case  with 
lens  is  0.95,  approximately  three  times  greater  than  without  the 
lens  (0.3),  which  is  an  indication  of  improved  sensitivity.  We 
move  the  target  along  the  lateral  axis  and  keep  the  antenna  and 


Fig.  6.  Maximum  differences  of  backscattered  signals  with  and  without  the 
lens.  The  target  is  a  metallic  object  (2.5  x  1.5  x  4.5  cm3).  At  each  lateral 
distance,  the  data  is  obtained  by  subtracting  the  two  peaks  in  the  difference  be¬ 
tween  the  backscattered  signal  with/without  the  presence  of  the  target  (Fig.  5). 
The  focus  is  at  zero  distance  from  the  planar  lens  surface. 


lens  static,  which  is  equivalent  to  shifting  the  imaging  system 
and  scanning  the  region  laterally. 

The  maximum  differences  for  the  cases  with  and  without  the 
lens  are  plotted  in  Fig.  6.  For  the  case  without  the  lens,  the 
maximum  differences  beyond  the  lateral  distances  of  ±6  cm 
are  in  actuality  the  same  as  the  noise  level,  which  makes  it  hard 
to  detect  the  target.  For  the  case  with  the  lens,  the  detectable 
range  is  approximately  the  same.  When  the  baseline  is  chosen 
as  0.15,  the  full  width  at  half-maximum  (FWHM)  for  the  case 
with  the  lens  is  approximately  4  cm,  whereas  the  FWHM  is  ap¬ 
proximately  6  cm  without  the  lens.  Since  the  FWHM  of  the 
scattering  difference  determines  the  image  resolution,  the  res¬ 
olution  is  therefore  increased  by  a  factor  of  1.5.  It  is  noticed 
that  although  the  focal  spot  size  is  reduced  by  a  factor  of  3  [7], 
the  resolution,  i.e.,  the  FWHM  of  maximum  difference,  is  not 
increased  accordingly  by  the  same  factor.  This  result  suggests 
that  the  focal  spot  size  is  not  the  only  factor  that  determines  the 
imaging  resolution,  as  the  receiving  characteristics  of  the  an¬ 
tenna  with  or  without  the  lens  also  contributes  to  the  difference 
in  the  received  signals. 

To  further  verify  the  imaging  resolution  of  the  lens-antenna 
system,  we  have  imaged  two  targets  of  equal  dimensions  (2.5  x 

1.5  x  4.5  cm3)  with  different  spacing.  Again,  the  two  objects  are 
shifted  from  one  edge  of  the  lens  to  the  other,  which  is  equivalent 
to  laterally  scanning  the  two  objects  with  the  antenna- with-lens 
system.  The  background  signals  are  subtracted  from  the  sig¬ 
nals  obtained  at  each  target  location.  The  value  of  the  difference 
is  used  to  construct  a  2-D  image  by  assigning  pixel  values  ac¬ 
cording  to  the  magnitude  of  the  signals  received.  Fig.  7  shows 
the  image  constructed  with  two  objects  with  spacing  of  3.5  and 

5.5  cm.  The  two  targets  are  irresolvable  in  the  lateral  direction 
when  the  spacing  is  shorter  than  3.5  cm.  However,  the  imaging 
system  is  able  to  distinguish  targets  for  a  spacing  of  5.5  cm,  so 
the  resolution  can  be  estimated  as  4  cm. 

While  we  have  demonstrated  confocal  scanning  in  the  lateral 
direction,  longitudinal  scanning  is  beneficial  to  probe  a  more 
deeply  buried  target.  Because  the  focus  is  fixed  along  the  longi¬ 
tudinal  axis  in  our  system,  the  confocal  configuration  can  only 
detect  targets  near  the  focus,  which  means  close  to  the  lens 
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Fig.  7.  Images  of  two  identical  metallic  objects  (2.5  x  1.5  x  4.5  cm3)  that 
are  separated  by  (a)  3.5  cm  and  (b)  5.5  cm.  The  two  objects  can  be  resolved  at 
a  lateral  spacing  of  5.5  cm,  but  are  irresolvable  for  3.5  cm. 


Longitudinal  distance  (cm) 

Fig.  8.  Maximum  differences  of  backscattered  signals  with  and  without  the 
lens.  The  target  is  a  metallic  object  ( 2 . 5  x  1 . 5  x  4 . 5  cm3 ) .  At  each  longitudinal 
distance,  the  data  is  obtained  by  subtracting  the  two  peaks  in  the  difference  be¬ 
tween  the  backscattered  signal  with/without  the  presence  of  the  target  (similar 
to  the  signals  shown  in  Fig.  5).  The  focus  is  at  zero  distance  from  the  planar  lens 
surface. 


surface.  Longitudinal  confocal  scanning  therefore  requires  dif¬ 
ferent  lenses  with  different  focus.  Here,  we  demonstrate  con- 
focal  detection  of  a  target  when  the  target  is  placed  along  the 
longitudinal  axis  of  the  lens.  Fig.  8  shows  the  maximum  dif¬ 
ferences  of  the  signals  in  the  cases  with  and  without  the  lens. 
At  the  focal  point  (longitudinal  distance  =  0),  the  difference  is 
greatest  and  is  about  three  times  higher  than  without  the  lens. 
Until  the  target  is  moved  2  cm  away  from  the  focus,  we  still 
obtain  a  higher  difference  than  without  the  lens.  Beyond  2  cm, 
the  difference  decreases  faster  with  the  lens  than  without,  as  the 


wave  exiting  the  lens  is  an  evanescent  wave  and  decreases  ex¬ 
ponentially  in  space. 


IV.  Conclusion 

We  have  explored  near-field  target  detection  with  a  focusing 
prolate  spheroidal  antenna  using  a  confocal  technique.  As  the 
antenna  was  designed  to  focus  in  free  space,  the  use  of  such  a 
system  is  limited  to  the  detection  of  targets  in  free  space  or  oth¬ 
erwise  embedded  in  a  weakly  scattering  medium.  When  linear 
targets  are  placed  at  the  focal  position,  the  antenna  is  sensitive 
to  their  orientation.  A  focusing  antenna  in  conjunction  with  a 
dielectric  lens  has  been  shown  to  have  a  resolution  of  approx¬ 
imately  4  cm.  The  resolution  can  be  further  improved  by  in¬ 
creasing  the  dielectric  constant  of  the  material  in  the  innermost 
layer  of  the  lens.  Our  present  confocal  system  prevents  a  scan¬ 
ning  of  depth  as  the  target  needs  to  be  placed  close  to  the  lens, 
within  2  cm  of  the  focus,  in  order  to  show  a  stronger  signal  re¬ 
turn  compared  to  without  the  lens.  Scanning  of  depth  may  be 
possible  if  the  probed  object  is  placed  within  a  lens  of  a  dif¬ 
ferent  design. 
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Abstract — This  paper  briefly  reviews  the  design  of  a  pro- 
late-spheroidal  impulse-radiating  hyperband  antenna  system  to 
launch  and  focus  fast  (100  ps)  high-voltage  (>100  kV)  pulses  onto 
biological  targets.  Design  and  experiments  on  a  graded  five-layer 
dielectric  lens,  designed  to  match  the  pulses  into  the  target 
medium,  are  outlined.  The  design  and  numerical  simulations  of 
a  switch  system  to  launch  high-voltage  spherical  TEM  pulses  are 
also  presented. 

Index  Terms — Apoptosis,  focusing  lens,  hyperband  antenna,  im¬ 
pulse-radiating  antenna,  launching  lens,  picosecond  pulses,  pro- 
late-spheroidal  reflector,  switch. 


I.  Introduction 

THE  EFFECT  of  high-intensity  nanosecond  electrical 
pulses  on  biological  cells  and  tissue  has  received  con¬ 
siderable  attention  by  the  scientific  community,  particularly 
over  the  past  five  years,  and  has  led  to  the  establishment  of  a 
new  research  field:  bioelectrics  [1].  In  the  time  domain  below 
approximately  100  ns,  the  rise  time  of  the  pulse  is  faster  than 
the  charging  time  of  the  plasma  membrane  of  most  mammalian 
cells,  meaning  that  the  field  will  pass  through  the  membrane 
into  the  cytoplasm.  High-voltage  pulses  in  the  nanosecond 
range  have  been  shown  to  penetrate  into  living  cells  to  per- 
meabilize  intracellular  organelles,  and  release  Ca2+  from  the 
endoplasmic  reticulum.  They  provide  a  new  approach  for 
physically  targeting  intracellular  organelles  with  many  appli¬ 
cations  including  activation  of  platelets  and  release  of  growth 
factors  for  accelerated  wound  healing  [2]  and  precise  control 
of  programmed  cell  death  (apoptosis),  which  has  been  shown 
to  cause  complete  elimination  of  melanoma  tumors  [3].  The 
pulses  were,  in  this  case,  delivered  by  needle-electrodes  to  the 
tumors. 
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The  bioelectric  effects  of  ultrashort  square  wave  pulses  were 
shown  to  follow  a  scaling  law  [4] 

S  =  S(EtN1^2)  (1) 

with  S  being  a  quantity  that  describes  the  intensity  of  the  observ¬ 
able  bioelectric  effect,  E  being  the  electric  field  amplitude,  r 
being  the  pulse  duration,  and  N  being  the  number  of  pulses.  This 
law,  which  was  shown  to  hold  for  nanosecond  pulses,  indicates 
that  for  picosecond  pulses  to  achieve  similar  effects  the  electric 
field  will  need  to  be  much  higher  than  for  nanosecond  pulses. 
However,  there  are  two  reasons  to  enter  the  picosecond-pulsed 
electric  field  range.  The  first  one  is  that  for  extremely  short 
pulses,  the  permittivity  of  the  various  cell  components,  rather 
than  their  resistivity,  determines  the  electric  field  distribution  in 
the  cell  [5].  The  electric  field  then  acts  directly  on  membrane 
proteins,  rather  than  causing  charging  of  the  membrane,  and, 
if  sufficiently  strong,  can  cause  direct  and  instant  conforma¬ 
tional  changes.  Subnanosecond  pulses  (200  ps)  were  found  to 
alter  the  cell  membrane  conductance  and  unirectifying  channels 
are  formed  when  cells  are  exposed  to  electric  field  strength  on 
the  order  of  20  kV/cm  [6]  for  2000  pulses.  The  disruption  of 
the  membrane  integrity  may  lead  to  the  change  of  physiological 
conditions  of  the  cell  and  cause  cell  death. 

Second,  besides  providing  the  opportunity  to  explore  a  new 
field  of  electric  field-cell  interactions,  subnanosecond  pulses 
will  ultimately  allow  medical  applications  for  delivery  of  pulsed 
electric  fields  without  invasive  electrodes,  using  antennas  in¬ 
stead.  A  possible  configuration  that  allows  us  to  generate  very 
high  electric  fields  uses  a  focusing  antenna  [7],  [8]  and  may  be 
able  to  induce  apoptosis  in  tissue  without  using  needles  as  the 
delivery  system  for  electric  pulses,  Fig.  1.  In  addition  to  using  a 
reflecting  antenna,  efforts  are  underway  to  utilize  lenses  in  com¬ 
bination  with  the  reflector  in  order  to  achieve  higher  spatial  res¬ 
olution  [9],  [10]. 

Besides  using  ultrashort  (10  ns)  electrical  pulses  for  medical 
therapies  [1 1],  the  pulsed  radiation  may  lead  to  medical  imaging 
methods.  Imaging  is  based  on  the  measurement  of  changes  in 
the  complex  permittivity  of  tissue,  and  may  complement  other 
methods  based  on  the  measurement  of  other  physical  parame¬ 
ters,  such  as  X-ray  computed  tomography  (CT),  magnetic  reso¬ 
nance  imaging  (MRI),  and  ultrasound.  Generally,  the  resolution 
is  limited  by  the  diffraction  time,  which  restricts  the  detectable 
target  size  to  larger  than  the  pulsewidth.  For  example,  the  spatial 
resolution  in  the  centimeter  range  is  expected  if  the  pulse  is  near 
100  ps.  In  such  an  imaging  system,  the  electromagnetic  waves 
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Fig.  1 .  Concept  of  treatment  of  cancerous  tissue  inside  the  body.  For  imaging  a 
second,  an  identical  IRA  in  a  confocal  arrangement  could  be  used  as  a  receiver. 


are  focused  on  a  target  inside  the  body.  The  scattered  signal  con¬ 
tains  information  on  the  dielectric  properties  and  the  geometry 
of  the  target  and  allows  its  identification  through  an  inverse  scat¬ 
tering  method.  3-D  scanning  allows  us  to  obtain  the  dielectric 
profile  of  the  tissue  and  to  detect  any  abnormality  in  a  uniform 
background. 

Research  on  the  use  of  impulse  radiating  antennas  (IRAs) 
for  medical  applications  requires,  in  addition  to  the  develop¬ 
ment  of  focusing  antennas,  research  on  the  biological  effect 
of  subnanosecond  pulses.  Experimental  studies  with  800-ps- 
long  pulses  with  electric  field  amplitudes  of  up  to  1  MV/cm  ap¬ 
plied  to  melanoma  (B16)  cells  have  provided  information  on 
the  pulse  parameters  required  to  induce  cell  death  [5].  Studies 
with  200-ps-long  pulses  are  underway.  First  results  indicate  that 
cell  death  requires  a  large  number  of  pulses  (>  106)  at  moderate 
electric  field  intensities  of  20  kV/cm  [6] .  However,  synergistic 
effects  such  as  caused  by  local  tissue  heating  might  allow  us  to 
reduce  the  pulse  number  considerably,  and  to  establish  IRAs  as 
a  new  tool  in  treatment  of  cancer. 

This  paper  presents  a  prolate- spheroidal  impulse  radiating 
antenna  (PSIRA)  system  to  focus  fast  (100  ps)  high-voltage 
(>100-kV)  pulses  on  to  a  biological  target.  The  treatment  of 
skin  cancer  (melanoma)  is  used  as  an  illustrative  example.  Cur¬ 
rently,  needle  arrays  are  used  for  treating  melanoma  tumors. 
This  approach  is  invasive  and  often  results  in  discomfort  to  the 
patient.  The  PSIRA  has  the  advantage  of  noninvasively  deliv¬ 
ering  the  required  pulsed  electric  fields  to  melanoma  tissues. 
By  directly  focusing  on  the  target  cells,  it  reduces  damage  to 
the  tissue  layers  surrounding  the  target  and  skin. 

The  design  of  the  PSIRA  and  focusing  lens  system  are  briefly 
reviewed  in  this  paper.  New  experimental  results  comparing  the 
focal  impulse  waveforms  and  the  beam  widths,  in  air  and  in  the 
focusing  lens,  are  presented.  It  is  found  that  the  focusing  lens 
can  be  viewed  as  a  bandpass  filter.  The  design  of  a  switch  system 
to  launch  high-voltage  (>100-kV)  spherical  TEM  pulses  from 
such  an  antenna  are  detailed.  The  experimental  results  on  the  fo¬ 
cusing  lens  and  the  design  of  the  switch  system  are  significant 
novel  contributions  that  complete  the  conceptual  realization  of 
the  PSIRA  system.  The  details  provided  in  this  paper  can  be 


used  for  the  fabrication  and  deployment  of  a  prototype  PSIRA 
system,  with  the  focusing  and  launching  lenses,  as  a  less  inva¬ 
sive,  safer,  and  more  effective  treatment  method. 

II.  PSIRA 

The  IRA  is  a  focused  aperture  hyperband1  antenna  suited  for 
radiating  very  fast  high-voltage  (>100-kV)  pulses  in  a  narrow 
beam  [14].  A  fast-rising  step-like  input  into  the  antenna  gives 
an  approximate  delta-function  response.  IRAs  thus  provide  an 
attractive  means  to  deliver  electric  pulses  to  induce  apoptosis  in 
tissue  cells  without  using  needles. 

A.  Prolate -Spheroidal  Reflector 

IRAs  are  composed  of  two  main  parts:  a  conical  TEM  trans¬ 
mission  line,  called  the  feed  arms,  and  a  focusing  optic,  which  is 
usually  a  reflector  or  a  lens.  Paraboloidal  reflector  IRAs,  where 
the  radiated  impulse  is  focused  at  infinity,  are  by  far  the  most 
common,  as  they  have  been  used  in  a  variety  of  applications  such 
as  transient  radars,  buried  target  identication,  etc.  However,  for 
use  in  the  biological  application  described  above,  where  electric 
fields  need  to  be  focused  onto  a  target  in  the  near  field,  a  pro¬ 
late- spheroidal  reflector  must  be  used  [7],  [15]. 

The  front  and  side  views  of  a  four  feed-arm  PSIRA  are  shown 
in  Fig.  2.  The  feed-arms  are  oriented  at  60°,  as  this  amplifies  the 
field  intensity  by  a  factor  of  1 .606  compared  to  the  two-arm  con¬ 
figuration  [16].  The  dimensions  of  the  feed  arms  are  determined 
by  their  200- Q  pulse  impedances. 

The  PSIRA  system  has  two  foci.  An  inhomogeneous  spher¬ 
ical  TEM  wave  launched  on  guiding  conical  conductors  from 
one  focus  (— ^o)  is  converted,  by  a  double  stereographic  trans¬ 
formation,  to  a  second  (reflected)  inhomogeneous  spherical 
TEM  wave  propagating  toward  the  second  focus  (zo)  [17]. 

B.  Focal  Waveform 

For  an  ideal  step  input,  the  electric  field  waveform  at  the 
second  focal  point  is  shown  in  Fig.  3.  The  waveform  has  three 
main  parts:  the  pre-pulse,  impulse,  and  post-pulse.  The  pre¬ 
pulse  is  the  direct  radiation  from  the  feed.  It  has  a  low  mag¬ 
nitude  and  lasts  for  a  long  duration.  The  impulse  is  due  to  the 
fields  from  the  reflector  converging  at  the  focal  point.  It  lasts 
for  a  short  time  and  has  a  large  amplitude.  The  post-pulse  is 
mostly  due  to  diffraction  effects  and  the  fields  outside  the  aper¬ 
ture.  The  pre-pulse  and  impulse  dominate  the  early-time  (high 
frequency)  behavior  of  the  IRAs,  while  the  post-pulse  provides 
information  on  the  late-time  (low  frequency)  behavior.  The  pro¬ 
late- spheroidal  reflector  differentiates  the  input;  the  derivative 
of  the  step  input  thus  results  in  the  impulsive  part  of  the  wave¬ 
form.  It  is  the  impulse  that  induces  electroporosis  or  apoptosis 
in  the  cells,  and  therefore,  must  be  maximized. 

C.  Optimization  of  Reflector  Dimensions 

The  pre-pulse  reduces  the  amplitude  of  the  impulse.  It  is 
therefore  desirable  to  have  as  small  a  pre-pulse  as  possible.  The 
magnitude  of  the  pre-pulse  is  solely  determined  by  the  size  of 
the  reflector  [15],  [17].  The  reflector  geometry  and  position  of 

Hhe  term  hyperband  refers  to  a  band  ratio  (not  bandwidth)  greater  than  one 
decade  [12].  It  is  a  well-established  International  Electrotechnical  Commission 
(IEC)  standard  (61000-2-13)  [13]. 
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After  investigation  of  various  electromagnetic  parameters, 
the  following  dimensions  were  found  to  be  most  suitable  [7], 
[15]: 

zp  =  0  a  =  0.625  m  6  =  4>  =  0.5  m  z0  =  0.875  m  /  =  !  m 

(2) 

where 

Zp  —  z  —  coordinate  of  the  truncation  plane 
a,  b  =  radii  of  the  prolate- spheroid 
zq  =  ( a 2  —  62)1/2  focal  distance 
l  =  a  +  zo  =  distance  used  for  normalization 
^  =  radial  coordinate 

as  shown  in  Fig.  2(a).  For  the  values  in  (2),  the  reflector  can  be 
easily  fabricated.  The  symmetry  of  the  geometry  greatly  also 
simplifies  the  problem. 

A  pulsewidth  of  100  ps  is  used  to  validate  the  applicability  of 
the  PSIRA  to  deliver  pulses  in  the  picosecond  regime.  A  100-ps 
pulse  duration  leads  to  a  reflector  with  practically  acceptable  di¬ 
mensions.  For  pulses  less  than  100  ps,  the  reflector  can  be  made 
smaller.  For  pulsewidths  greater  than  100  ps,  a  larger  system 
is  required.  However,  using  the  PSIRA  system  for  pulses  with 
durations  much  greater  than  100  ps,  e.g.,  >500  ps,  is  not  recom¬ 
mended,  as  the  size  of  the  system  may  not  be  very  practical.  The 
reflector  dimensions  in  (2)  were  optimized  for  a  100-ps  ramp 
rising  input. 

CST  Microwave  Studio,  a  3-D  finite  integral  time  domain 
(FITD)  commercially  available  software,  was  used  for  simula¬ 
tions.  A  1-V  100-ps  ramp  rising  step  was  used  as  the  input  in  all 
simulations. 


Fig.  2.  Side  and  front  views  of  the  prolate- spheroidal  reflector,  (a)  Side  view, 
(b)  Front  view. 


the  truncation  plane  zp  are  optimized  to  obtain  an  acceptable 
compromise  between  the  physical  dimensions  of  the  reflector 
and  the  pre-pulse  amplitude. 


III.  Focusing  Lens 

Section  II  considered  the  design  of  a  prolate- spheroidal  re¬ 
flector  for  concentrating  a  pulse,  in  the  near  field,  from  one  of 
the  foci  onto  a  target  at  the  second  focus.  The  idea  is  to  get  a 
very  fast  very  intense  electromagnetic  pulse  to  illuminate  the 
target  (e.g.,  a  tumor).  By  very  fast  we  mean  fast  enough  to  get 
the  spot  diameter  (beamwidth)  at  the  focus  down  to  the  target 
size,  or  as  close  to  this  as  technology  allows.  At  the  same  time, 
a  small  spot  size  also  implies  large  fields  [18]. 

One  problem  with  placing  fields  on  the  target  concerns  the 
dielectric  properties  of  the  target  medium  and  its  surroundings. 
If  the  wave  incident  on  the  target  is  in  air,  but  the  target  medium 
has  a  large  relative  permittivity  (say,  about  81  for  water),  then 
there  will  be  a  significant  reflection  of  the  pulse,  leading  to  a 
smaller  field  in  the  target  medium  [18].  This  section  discusses 
the  addition  of  a  lens,  called  the  focusing  lens,  to  better  match 
the  wave  to  the  target,  and  to  reduce  spot  size,  thereby  increasing 
the  field  on  the  target.  A  schematic  of  the  PSIRA  and  focusing 
lens  system  is  shown  in  Fig.  4.  A  brief  overview  of  the  design 
and  experiments  on  a  five-layer  lens  is  presented  [10]. 

The  PSIRA  and  lens  system  can  deliver  100-ps  electric  field 
pulses  with  a  spatial  resolution  in  the  millimeter  range.  While 
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Fig.  4.  PSIRA  and  lens  system  to  focus  and  match  into  (biological)  target 
medium  [10]. 


our  experiments  use  10-V  45-ps  rising  signals,  the  design  is  for 
100  kV  or  greater  to  achieve  the  desired  field  strengths  with 
100-ps  rising  pulses. 

A.  Summary  of  Design 

The  design  of  a  spatially  limited  graded  hemispherical  di¬ 
electric  focusing  lens,  with  relative  permitivity  varying  from 
ermin  =  1  to  ermax  =  81,  is  based  on  a  transmission  line  model. 
There  are  two  important  parameters  to  consider  in  the  design, 
which  are:  1)  the  dimension  of  the  lens  and  2)  the  number  of 
layers.  Detailed  calculations  to  minimize  the  lens  radius  and  op¬ 
timize  the  number  of  layers  are  provided  in  [9]  and  [10].  The 
lens  is  designed  such  that  all  layers  have  the  same  electrical 
width.  A  ten-layer  lens  with  a  15 -cm  radius  was  found  to  pro¬ 
vide  satisfactory  electric  enhancement. 

Based  on  plane- wave  approximations,  scaling  relationships 
for  various  electromagnetic  parameters  are  calculated  in  [10] 
and  [19].  The  following  results  are  most  significant: 

•  electric  enhancment  =  e^ax ; 

•  focal  spot  size  reduction  =  e^/x2 ; 

•  displacement  enhancement  =  Lax ; 

as  they  are  used  to  quantify  the  efficacy  of  the  lens.  The  large 
displacement  enhancement  is  extremely  advantageous,  as  it  is 
the  electric  displacement  that  acts  on  the  (abnormal)  cells. 

The  graded  lens  design  assumes  that  the  biological  target  is 
placed  in  direct  contact  with  the  last  layer.  Some  preliminary 
considerations  regarding  the  problems  and  design  options  for 
focusing  more  deeply  buried  targets  are  presented  in  [20] .  The 
focal  spot  size  will  be  affected  if  the  lens  is  used  on  an  uneven 
surface  as  the  rays  arrive  at  the  focal  point  at  different  times.  For 
such  uneven  surfaces,  more  innovative  techniques  for  using  the 
lens  are  currently  being  explored  [10].  As  a  proof-of-concept,  a 
five-layer  focusing  lens,  er  .  =  1  to  er  =  9,  is  considered 
in  this  paper.  Numerical  simulations  and  experiments,  explained 
in  Section  IV,  are  used  to  demonstrate  the  viability  and  applica¬ 
bility  of  such  a  system. 


feed  arms 


mnJfxJronix  TDS  8000B 

Oscilloscop, 


prolate-spheroidal 

reflector 


— 

5-layer  hemispherical 
focusing  lens 


9 


D-aot  probe 


Picosecond  Pulse  Labs 
Pulse  generator 
45  ps  10  V 


Fig.  5.  Experimental  setup  of  60°  four  feed-arm  PSIRA  with  focusing  lens, 
target,  and  D-dot  probe  [10]. 


pulse  generator 

I 


pulse  sharpener  1 


PSIRA 

(first  focal  point) 


Focusing  lens 


Target  medium  + 
D-dot  probe 
(second  focal  point) 


[  100  ps  filter  ] 


sampler  1 


sampling 
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Fig.  6.  Block  diagram  of  the  experimental  setup. 


B.  Experimental  Results 

Fig.  5  shows  the  experimental  setup.  A  block  diagram  of  the 
setup  is  shown  in  Fig.  6.  It  consists  of  a  prolate- spheroidal  re¬ 
flector  and  two  feed  arms  at  60°  over  a  ground  plane,  i.e.,  a  60° 
four  feed-arm  PSIRA  configuration.  A  10-V  45-ps  ramp-rising 
signal,  generated  by  a  Picosecond  Pulse  Laboratories,  Boulder, 
CO2  pulse  generator,  was  used  as  a  source  at  the  first  focal  point. 
Our  analytical  calculations  assume  an  input  with  a  rise  time  of 
100  ps.  Therefore,  a  filter  is  placed  between  the  D-dot  probe 
and  the  sampler.  This  filter  converts  the  45-ps  signal  from  the 
pulser  to  one  with  a  rise  time  of  100  ps  [21].  The  lens  and  target 
medium  (slab)  are  placed  at  the  second  focal  point. 

The  five-layer  hemispherical  focusing  lens  and  target  were 
made  to  our  specifications  by  TPL  Inc.,  Albuquerque,  NM.3 
The  nanocomposite  dielectric  material  developed  by  TPL  Inc. 
is  ideal  for  the  fabrication  of  our  multilayer  lens  since  desired 
dielectric  constants  can  be  achieved.  The  lens  is  placed  at  the 
second  focal  point.  A  slab  placed  in  front  of  the  lens  is  used 
to  simulate  the  target  medium.  A  handmade  fast  D-dot  sensor 
is  placed  inside  the  target  medium.  The  signal  from  the  D-dot 
probe  was  measured  using  a  Textronix  TDS  8000B  digital  sam¬ 
pling  oscilloscope.  The  experimental  results  are  normalized  to 
a  1-V  input  [10]. 

Complete  details  of  the  D-dot  probe  are  provided  in  [10]  and 
[21].  Basically,  the  probe  is  a  3-mm-high  2-D  teardrop-shaped 
electrode  designed  to  have  a  pulse  impedance  of  50  D  in  the 
dielectric  medium,  i.e.,  150  O  in  air.  The  impedance  mismatch 

2[Online].  Available:  http://www.picosecond.com/ 

3 [Online].  Available:  http://www.tplinc.com/ 
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Fig.  7.  Comparison  of  focal  impulse  waveforms,  at  the  second  focal  point,  in 
air  and  inside  the  lens,  with  the  100-ps  filter. 


Distance  from  focal  point  (cm) 

Fig.  8.  Comparison  of  simulation  and  experiment  results  of  beamwidth  in  air 
and  inside  the  lens. 


between  the  probe  tip  and  the  input  pulse,  for  the  measurements 
in  air,  causes  reflections.  However,  the  probe,  of  height  h  = 
3  mm,  is  resonant  at  frequencies  above 


/ 


c 

n 

4 


3  x  1011  mm 
0.75  mm 


0.4  THz 


(3) 


i.e.,  in  a  quasi-static  sense,  the  measurements  are  valid  over  a 
very  large  frequency  spectrum  of  the  input  pulse  and  are  there¬ 
fore  considered  to  be  sufficiently  accurate. 

The  measurement  of  fast  (<  100  ps)  electromagnetic  pulses 
inside  a  high-permittivity  medium  presents  some  challenges. 
These  are  especially  associated  with  the  shorter  wavelengths 
and  the  contact  of  the  D-dot  antenna  with  the  medium.  Some 
techniques  to  mitigate  these  problems  to  some  degree  are  dis¬ 
cussed  extensively  in  [22]-[24]. 

Fig.  7  compares  the  focal  impulse  waveforms  in  air  and  in¬ 
side  the  focusing  lens.  The  amplitude  of  the  impulse  in  air  is 
5.44  V/m  and  the  full  width  at  half  maximum  (FWHM)  is  85  ps. 
One  would  expect  the  FWHM  to  be  identical  to  the  maximum 
rate  of  rise  of  the  input  pulse,  i.e.,  tmr  =  100  ps  [21].  The 
reason  for  this  discrepancy  is  that  the  higher  frequencies  in  the 
input  pulse  are  better  focused  than  the  lower  frequencies,  i.e., 
the  focusing  action  itself  filters  out  lower  frequency  content  in 
the  focused  waveform.  Similar  results  observed  in  simulations 
corroborate  this  explanation.  The  amplitude  of  the  impulse  in 
the  focusing  lens  is  7.66  V/m  and  the  FWHM  is  127.5  ps.  The 
large  FWHM  is  most  likely  due  to  loss  and  dispersion  in  the 
lens  materials  (without  loss  and  dispersion,  the  FWHM  would 
be  100  ps). 

In  Fig.  7,  one  observes  the  enhancement  inside  the  lens.  Since 
the  half-widths  of  the  impulse  responses  in  air  and  in  the  lens  are 
not  identical,  the  electric  enhancement  is  given  by  the  ratio  of 
the  areas,  A  =  (halfwidth)  (peak  amplitude)  =  FWHM-Fmax, 
under  the  curves,  i.e., 


^-impulse-air  —  ($5  ps)  (5.44  V/m) 

=  462.4  ps  V/m 

^impulse-lens  =  (127.5  ps)  (7.66  V/m) 
=  976.65  ps  V/m. 


TABLE  I 

Beamwidth  of  Electric  Field  in  Air  and  With  Focusing  Lens 


Simulation 

Experiment 

No  lens  (EFSSnl  in  cm) 

3.61 

4.451 

With  Lens  (EFSSwl  in  cm) 

1.19 

2.036 

EFSSwl /EFSSnl 

3.04 

2.187 

The  electric  enhancement  is  impulse- lens Mimpuise- air  ~  2.11. 
This  is  much  larger  than  the  analytical  estimate,  ej4  — 
8.21/4  «  1.69.  The  reason  for  this  disagreement  is  due  to  the 
85 -ps  FWHM  of  the  impulse  in  air  [21].  The  analytical  approx¬ 
imations  do  not  consider  the  filtering  of  the  lower  frequencies 
due  to  the  focusing  of  the  wave,  at  the  second  focal  point,  in 
air.  This  results  in  a  much  higher  electric  enhancement,  as 
observed  above.  Of  course,  one  could  also  simply  examine  the 
ratio  of  the  peak  electric  fields  inside  the  lens  and  in  air,  which 
is  7.66/5.44  =  1.408.  The  analogous  enhancement  obtained 
from  simulations  is  10.725/6.247  =  1.717. 

Fig.  8  compares  the  numerical  simulation  and  experimental 
results  of  the  beamwidth  in  air  and  in  the  lens.  Errors  in  the 
experimental  data  are  of  the  order  of  the  spot  size.  Note  that 
the  geometric  focal  point  in  the  lens  is  shifted  to  the  right  by 
approximately  x  =  +0.25  cm.  The  discrepancy  between  the 
experimental  and  numerical  results  is  due  to  losses  occurring 
in  the  lens  materials;  the  width  of  the  pulse  emerging  from  the 
lens  is  127.5  ps,  which  results  in  a  larger  spot  size.  A  summary 
of  the  spot  sizes  is  given  in  Table  I.  Again,  these  results  agree 
well  with  those  expected  from  the  analytical  calculations. 

C.  Focusing  Lens  as  a  Bandpass  Filter 

The  ratio  of  the  magnitudes  of  the  Fourier  transform  of  the 
impulse  waveforms  from  experiments,  in  the  lens  and  in  air,  is 
shown  in  Fig.  9. 

The  bandwidth  of  the  lens  is  approximately  3  GHz,  from  0.4 
to  3.4  GHz.  Above  3.5  GHz,  there  is  almost  no  electric  field  en¬ 
hancement  (analytical  or  experimental).  The  focusing  lens  thus 
acts  as  a  bandpass  filter,  filtering  frequencies  below  0.4  GHz 
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Frequency  (GHz) 

Fig.  9.  Ratio  of  magnitude  of  Fourier  transform  of  the  focal  impulse  waveform 
from  experiments. 


Fig.  10.  Ratio  of  magnitude  of  Fourier  transform  of  the  focal  impulse  wave¬ 
form,  from  simulations,  for  an  input  pulse  with  100-ps  rise  time.  The  dotted  gray 
line  is  1.77,  the  electric  enhancement  obtained  with  the  lens. 


and  above  3.4  GHz.  The  filtering  of  higher  frequencies  is  due  to 
losses  occurring  in  the  lens  materials,  while  the  the  lower  fre¬ 
quencies  are  filtered  due  to  the  lens  geometry  [25]. 

The  curve  in  Fig.  9  is  compared  to  similar  results  obtained 
from  simulations  in  Fig.  10.  The  (average)  amplification  is 
approximately  1.77  (dotted  gray  line)  up  to  a  frequency  of 
10  GHz,  as  expected,  since  the  lens  materials  are  assumed 
lossless  and  dispersionless  in  the  simulations.  This  is  closer  to 
that  obtained  from  the  ratio  of  the  peak  focal  electric  fields,  i.e., 
10.725/6.247  =  1.717.  Therefore,  the  analogous  experimental 
enhancement,  7.66/5.44  =  1.408,  is  considered  in  Fig.  9. 

IV.  Switch  System  Design 

For  input  voltages  of  100  kV  or  more,  a  switch  system  and 
“launching”  lens  are  necessary  to  effectively  launch  a  spherical 
TEM  wave  from  the  first  focal  point.  The  switch  system  consists 
of  switch  cones,  a  pressure  vessel  (PV),  and  a  gas  (typically 
hydrogen)  chamber.  The  hydrogen  chamber  (HC)  is  required 
to  support  the  high  field  intensities  without  breaking  down,  and 
hence,  allow  conduction.  The  design  and  numerical  simulations 
of  such  a  switch  system  and  launching  lens  are  presented  in 
this  section.  Although  the  input  for  the  numerical  simulations 


is  only  1  V,  analogous  designs  for  the  prototype  IRA  tested  at 
high- voltages  (>100  kV)  [26]  indicate  that  similar  results  can 
be  expected  for  the  switch  system  configurations  in  this  section. 

A.  Initial  Launching  Lens  Designs 

The  initial  approach  to  the  design  of  a  launching  lens  followed 
methodologies  similar  to  those  used  in  the  prototype  IRA  [27]. 
This  approach  assumes  the  feed  point  and  the  focal  point  to  be 
spatially  isolated.  The  objective  of  the  lens  design  is  to  ensure 
that,  within  the  lens,  a  spherical  TEM  wave  is  centered  on  the 
feed  point.  However,  outside  the  lens,  an  approximate  spher¬ 
ical  TEM  wave  is  centered  at  the  first  focal  point.  Uniform  and 
nonuniform  designs  explored  are  summarized  below. 

1)  Uniform  Lens  Design:  Analytical  calculations  defining 
the  boundary  of  a  dielectric  lens,  with  fixed  relative  permittivity, 
are  detailed  in  [28].  These  equations  are  derived  using  high-fre¬ 
quency  (optical)  approximations.  It  is  shown  that  the  dielectric 
constant  required  for  such  a  design  must  be  greater  than  25.  This 
constraint  is  a  serious  drawback  since  materials  with  such  high 
dielectric  constants  have  a  large  loss  and  dispersion. 

2 )  Nonuniform  Lens  Design:  Lower  dielectric  constants  can 
be  obtained  by  designing  a  lens  in  which  the  er  varies  (dis¬ 
cretely)  across  the  lens  boundary.  Note  that  the  lens  is  assumed 
to  be  a  body-of-revolution.  An  analytically  aided  simulation  ap¬ 
proach  is  used.  The  lens  boundary,  dielectric  constants  of  var¬ 
ious  layers,  and  thicknesses  of  the  layers  are  first  estimated  an¬ 
alytically.  Simulations  are  then  used  to  examine  the  time-of-ar- 
rival  of  electric  fields  on  a  spherical  surface  to  ensure  that  a 
spherical  TEM  wave  originates  from  the  first  focal  point.  This 
procedure  is  iterated  until  the  design  is  achieved  within  the  de¬ 
sired  tolerance.  Two  designs  are  described  below:  the  planar  and 
the  conical  [29]-[32]. 

Planar  Design:  The  planar  lens  is  one  in  which  the  rela¬ 
tive  permittivity  of  the  layers  varies  along  the  (rotational)  axis 
of  the  lens.  The  design  is  log  periodic,  i.e.,  all  layers  have  the 
same  electrical  length.  Three-  and  six-layer  designs  were  inves¬ 
tigated.  Simulations  indicated  that  the  time  spread  in  the  rays 
originating  from  the  source  was  over  100  ps.  This  time  spread 
is  highly  undesirable,  as  it  indicates  that  the  spherical  wavefront 
is  distorted  well  beyond  the  acceptable  tolerance  of  20  ps. 

Conical  Design:  In  this  design,  the  relative  permittivity  of 
the  layers  is  a  function  of  the  polar  angle  8.  Simulation  results 
from  a  seven-layer  configuration  were  almost  identical  to  the 
planar  design. 

Due  to  the  very  large  distortion  of  the  spherical  wave,  the 
planar  and  conical  designs  were  not  iterated,  as  this  would  have 
been  too  resource  intensive.  The  large  distortion  in  the  spherical 
wavefront  could  be  due  to  many  reasons,  such  as  the  optical 
approximations  breaking  down  at  lower  frequencies,  reflections 
from  inner  layers,  etc.  This  approach  was  abandoned  in  favor  of 
the  simpler  designs  described  in  Sections  V  and  VI. 

B.  Investigation  of  Various  Switch  Configurations 

In  the  design  of  the  launching  lens  described  above  it  was  as¬ 
sumed  that  the  feed  point  and  first  focal  point  were  spatially  iso¬ 
lated.  This  assumption  imposed  a  severe  constraint  on  the  lens 
design.  The  problem  is  greatly  simplified  if  the  geometric  center 
of  the  feed  point  is  made  to  coincide  with  the  first  focal  point. 
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TABLE  II 

Peak  Focal  Impulse  Amplitude  and  Beamwidth  for 
Various  Switch  Configurations 


Configuration 

Fmax  (V/m) 

Spot  size  (cm) 

VBCS-75 

6.727 

3.894 

VBCS-100 

5.619 

4.065 

VBCS-125 

4.708 

4.304 

VBCS-150 

4.057 

4.414 

4FASC 

6.363 

3.546 

T4FASC 

7.217 

3.889 

S4CFA 

5.114 

3.39 

VBCS-Z  denotes  a  vertical  bicone  switch  with  impedance  2 Z  Q. 

Before  proceeding  to  the  launching  lens  designs,  various  switch 
configurations  are  described  where  the  geometric  center  of  the 
switch  cones  is  the  first  focal  point.  The  focal  impulse  ampli¬ 
tude,  Tmax?  and  beamwidth  were  compared  for  the  following 
switch  configurations  [33]. 

1)  Four  Feed  Arms  with  Switch  Cones  (4FASC):  A  200-Q  bi¬ 
cone  switch,  of  height  1 .0  cm,  centered  at  the  first  focal 
point  transitions  to  the  feed-arms  of  the  PSIRA.  The  con¬ 
nection  between  the  switch  cone  base  and  the  feed-arms  is 
called  the  loft  connection. 

2)  Truncated  Four  Feed  Arms  with  Switch  Cones  (T4FASC): 
This  configuration  is  identical  to  the  4FASC  except  that  the 
feed  arms  are  truncated  at  a  distance  of  19  cm  from  the  first 
focal  point. 

One  notes  that  the  rise  time  of  tg  =  100  ps  corresponds  to  a 
physical  distance  of  ctg  =  3  cm  (in  air).  Therefore,  a  switch 
of  appropriate  dimensions  can  also  serve  as  a  guiding  structure. 
The  following  two  such  configurations  were  studied. 

1)  Vertical  Bicone  Switch  ( VBCS '):  A  vertical  bicone  is  used 
as  the  source  where  the  height  of  each  cone  is  6.0  cm.  This 
structure  also  serves  to  guide  the  spherical  TEM  waves 
originating  from  the  feed  point.  The  impedance  of  the 
cones  was  varied. 

2)  Slanted  Four-Cone  Feed  Arms  (S4CFA):  This  structure  is 
almost  identical  to  the  T4FASC,  except  that  the  feed  arms 
are  replaced  by  200-Q  cones.  Each  cone  is  6.0  cm  in  height. 

Results  from  numerical  simulations  for  ^max  and  the  spot 
size,  at  the  second  focal  point,  for  the  various  configurations  are 
summarized  in  Table  II. 

One  finds  the  VBCS-75  and  T4FASC  designs  to  be  the  most 
promising.  Both  these  configurations  are  easy  to  fabricate,  yield 
a  relatively  high  peak  electric  field  and  an  acceptable  spot  size 
compared  to  the  other  configurations  in  Table  II.  Additionally,  it 
is  also  found  that  the  pre-pulse  is  shorter  with  these  structures, 
compared  to  those  in  [10].  Therefore,  these  configurations  are 
further  considered  for  the  launching  lens  designs  in  Section  V. 

C.  PV  Designs 

The  gas  chamber,  PV,  and  launching  lens  are  necessary  as, 
in  the  final  experimental  realization  of  the  PSIRA,  100  kV  or 
more  will  be  applied  across  the  switch  gap.  The  gas  chamber 
typically  contains  hydrogen  or  SFe  under  high  pressure.  In  this 
paper,  hydrogen  is  used  as  an  example.  As  mentioned  in  Section 
IV-B,  the  location  of  the  geometric  center  of  the  switch  cones 
at  the  first  focal  point  greatly  simplifies  the  design.  The  design 


Fig.  1 1 .  Side  view  of  the  T4FASC  configuration  with  a  spherical  PV  and  spher¬ 
ical  HC. 


TABLE  III 

Dimensions  of  Switch  System  Components  for 
the  T4FASC-SPV  Design 


Component 

Height  (cm) 

Radius  (cm) 

switch  cone 

h  =  0.5 

0.51 

cylindrical  support 

-Hess  —  0.5 

0.51 

hydrogen  chamber 

0.71 

pressure  vessel 

- 

7pv  =  2.0 

oil  medium 

— 

r0ii  =  5.0 

is  further  simplified  if  one  considers  the  PV  to  serve  the  func¬ 
tion  of  the  lens.  The  switch  system  then  consists  of  only  three 
components,  which  are:  1)  switch  cones  and  guiding  structures; 
2)  HC;  and  3)  PV. 

The  PV  may  be  cylindrical  or  spherical.  The  peak  focal  im¬ 
pulse  amplitudes  of  the  VBCS-75  configurations,  with  a  spher¬ 
ical  PV,  were  approximately  30%  lower  than  the  corresponding 
T4FASC  designs.  Therefore,  numerical  simulation  results  for 
only  the  T4FASC  configuration  with  the  spherical  and  cylin¬ 
drical  PV  designs  are  presented  here  [34],  [35]. 

1 )  T4FASC  With  a  Spherical  PV  ( T4FASC-SPV ):  To  provide 
structural  support  to  the  PV,  a  cylindrical  support  section  (CSS), 
of  height  Hcss,  is  added  to  the  T4FASC  design.  The  PV  and  HC 
are  spherical.  The  relative  permittivity  of  the  PV  medium  is  as¬ 
sumed  to  be  3.7.  The  impedance  of  the  switch  cones  is  200  Cl  in 
the  PV  medium,  i.e.,  half-angle  0  =  45.58°.  A  spherical  con¬ 
tainer,  with  er  =  2.25  (transformer  oil),  surrounds  the  PV  and 
is  used  to  denote  the  oil  “bath”  that  would  be  used  in  the  final 
practical  design.  Fig.  1 1  shows  the  side  view  of  the  various  com¬ 
ponents  of  the  switch  system.  Dimensions  of  these  components 
are  summarized  in  Table  III. 

Probes  to  monitor  the  time-of-arrival  of  electric  fields  were 
placed  on  a  (virtual)  sphere  of  radius  10  cm  (near  field).  The 
normalized  responses  from  these  probes  are  shown  in  Fig.  12 


KUMAR  et  al.\  HYPERBAND  ANTENNA  TO  LAUNCH  AND  FOCUS  FAST  HIGH-VOLTAGE  PULSES  ONTO  BIOLOGICAL  TARGETS 


1097 


(a) 


l 

. .  .  1  '  V 

<#'  ; 

E 

> 

0.8 

M  \ 

CD 

c 

CO 

0.6 

M  : 

CL 

£ 

0.4 

M  - 

«x> 

LD 

0.2 

n 

- 

i  .  ,—lrV^  r  .  i  .  .  1 . 1 . 1  ' 

350  400  450  500  550  600 

Time  (ps) 

0°  -  45°  90° 

30°  - . .  60° 

(b) 


(b) 


(c) 


Time  (ps) 

0°  -  45°  90° 

30°  -  60° 

(d) 

Fig.  12.  Normalized  Ee  and  E $  components  of  the  responses  from  the  electric 
field  probes  on  the  xy-,  —  yz~,  and  —zx -planes  for  the  T4FASC-SPV  configu¬ 
ration.  (a)  Normalized  E $  in  the  xy- plane  (b)  Normalized  Ee  in  the  xy- plaen 
(c)  Normalized  Ee  in  the  —yz- plane  (d)  Normalized  E $  in  the  —zx- -plane. 


(each  response  is  normalized  with  respect  to  its  maximum).  The 
component  in  the  —yz- plane  and  the  Eq  component  in  the 
—zx- plane  are  0.  The  time  spread  in  the  electric  fields  is  less 
than  20  ps,  which  is  within  the  acceptable  tolerance. 


Fig.  13.  Electric  field  focal  impulse  waveform  for  the  T4FASC-SPV  configu¬ 
ration. 


Distance  from  focal  point  (cm) 

Fig.  14.  Spot  size  for  the  T4FASC-SPV  configuration. 

The  focal  waveform  and  spot  size  are  shown  in  Figs.  13  and 
14.  These  results  indicate  that  the  electric  field  is  enhanced, 
compared  to  the  T4FASC  configuration  in  Table  III,  with 
the  PV  by  approximately  91%,  ^max  =  13.76  V/m,  with 
a  corresponding  increase  of  only  19%  in  the  spot  diameter, 
beam  width  =  4.64  cm.  The  electric  field  is  enhanced  due 
propagation  through  the  PV  dielectric  medium. 

2)  T4FASC  With  a  Cylindrical  PV  (T4FASC-CPV):  Analyt¬ 
ical  calculations  using  optical  approximations  for  a  T4FASC 
configuration  with  a  cylindrical  PV  and  cylindrical  HC  are  de¬ 
tailed  in  [36].  It  is  shown  that  for  a  given  relative  permittivity, 
radius,  and  height  of  the  PV,  the  surrounding  oil  medium  can  be 
used  as  the  launching  lens.  Fig.  15  shows  the  side  view  of  the 
various  components  of  the  switch  system. 

The  dimensions  of  the  components  are  summarized  in 
Table  IV.  The  radius  and  height  of  the  PV  are  determined  using 
the  formulas  in  [36]  for  ryc  =  ht&nO  =  0.5  tan(45.58°)  = 
0.51  cm  and  eril  =  2.25;  eril  is  the  relative  permittivity  of  the 
launching  lens.  Note  that  hpv  =  h  +  Hcss  in  Table  IV. 

The  time-of-arrival  of  the  electric  fields  on  probes  placed  on 
a  spherical  surface,  of  radius  10  cm,  showed  a  time  spread  sim¬ 
ilar  to  that  of  the  T4FASC-SPV,  i.e.,  less  than  20  ps.  The  focal 
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Fig.  15.  Side  view  of  the  T4FASC  configuration  with  a  cylindrical  PV  and 
cylindrical  HC. 


TABLE  IV 

Dimensions  of  Switch  System  Components  for 
the  T4FASC-CPV  Design 


Component 

Height  (cm) 

Radius  (cm) 

switch  cone 

h  =  0.5 

0.51 

cylindrical  support 

f7css  =  0.545 

0.51 

hydrogen  chamber 

Zihc  =  1-0 

0.51 

pressure  vessel 

V  =  1.045 

rpv  =  1.905 

oil  medium 

— 

Gfil  —  5.0 

Fig.  16.  Electric  field  focal  impulse  waveform  for  the  T4FASC-CPV  configu¬ 
ration. 

waveform  and  spot  size  are  shown  in  Figs.  16  and  17.  These 
results  indicate  that  the  electric  field  is  enhanced,  compared  to 
the  T4FASC  configuration  in  Table  IV,  with  the  PV  by  approxi¬ 
mately  62%,  Hmax  =  11.7  V/m.  The  corresponding  increase  in 
the  spot  diameter  is  only  18%,  beamwidth  =  4.95  cm.  Although 
Hmax  is  less  than  T4FASC-SPV,  the  T4FASC-CPV  configura¬ 
tion  is  attractive  since  it  is  easier  to  fabricate. 

Various  components  of  the  switch  and  PV  design  were  opti¬ 
mized  so  that  a  large  peak  electric  field,  with  a  small  spot  size, 
was  obtained.  For  the  spherical  PV  configuration,  the  following 
ranges  yield  a  reasonable  electric  field  amplification: 


Distance  from  focal  point  (cm) 


Fig.  17.  Spot  size  for  the  T4FASC-CPV  configuration. 

•  feed  arm  length:  18  cm  <  /  <  22  cm; 

•  height  of  cylindrical  support  structure:  0.2  cm  <  Hcss  < 
0.6  cm; 

•  PV  radius:  1.50  cm  <  rpv  <  3.0  cm. 

The  height  and  radius  of  the  PV  for  the  cylindrical  PV  design 
are  fixed  by  the  formulas  in  [36].  The  range  of  the  feed  arm 
length,  Z,  for  the  T4FASC-CPV  configuration  is  identical  to  the 
T4FASC-SPV  design. 

V.  Discussion 

The  design  of  the  PSIRA,  with  the  focusing  lens  and  switch 
system,  presented  above  is  conceptually  complete.  The  system 
can  deliver  <200-ps  pulses  with  a  spot  area  of  approximately 
2  cm2,  which  is  adequate  for  biological  applications.  However, 
there  remain  two  major  technological  challenges  before  a  pro¬ 
totype  device  can  be  field-tested. 

The  first  concerns  the  dielectric  materials  used  in  the  focusing 
lens.  As  observed  in  Section  III,  the  amplitude  of  the  impulse 
is  reduced,  and  the  beamwidth  increased,  due  to  losses  occur¬ 
ring  in  the  lens  materials.  It  is  also  observed  that  the  lens  filters 
frequencies  above  5  GHz.  Therefore,  materials  with  less  losses 
and  a  better  frequency  response,  from  0  to  10  GHz,  must  be  ex¬ 
plored. 

The  second  hurdle  pertains  to  high-frequency  switching.  Nu¬ 
merical  simulations  used  in  the  design  of  the  switch  system  in 
Section  IV  do  not  consider  the  physics  of  gas  discharges.  For  the 
switch  designs  considered,  current  technology  limits  the  repe¬ 
tition  rate  to  a  few  hundred  hertz.  Higher  frequencies,  of  the 
order  of  a  few  kilohertz,  are  required  to  induce  apoptosis  in 
cells.  Hence,  techniques  to  increase  the  frequency  of  switch  dis¬ 
charges  must  be  investigated. 

VI.  Concluding  Remarks 

This  paper  has  reviewed  the  design  of  a  prolate- spheroidal 
impulse-radiating  hyperband  antenna  to  concentrate  100-ps 
pulses  onto  biological  targets.  The  design  of  a  graded  five-layer 
focusing  lens  system,  to  better  match  the  pulses  from  air  into 
the  target  medium,  has  been  outlined.  Experimental  results 
with  the  PSIRA  and  lens  system  agree  well  with  those  obtained 
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from  numerical  simulations.  It  has  been  shown  that  the  fo¬ 
cusing  lens  acts  as  a  bandpass  filter.  The  design  and  numerical 
simulations  of  a  switch  system  have  been  presented.  Spherical 
(T4FASC-SPV)  and  cylindrical  (T4FASC-CPV)  PV  designs 
have  been  shown  to  yield  an  enhanced  focal  impulse  with  less 
than  20-ps  distortion  in  the  spherical  wave  originating  from  the 
first  focal  point. 
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Subnanosecond  Electric  Pulses  Cause  Membrane 
Permeabilization  and  Cell  Death 

Shu  Xiao*,  Siqi  Guo,  Vasyl  Nesin,  Richard  Heller,  and  Karl  H.  Schoenbach 


Abstract — Subnanosecond  electric  pulses  (200  ps)  at  electric  field 
intensities  on  the  order  of  20  kV/cm  cause  the  death  of  B16.F10 
murine  melanoma  cells  when  applied  for  minutes  with  a  pulse  rep¬ 
etition  rate  of  10  kHz.  The  lethal  effect  of  the  ultrashort  pulses  is 
found  to  be  caused  by  a  combination  of  thermal  effects  and  electri¬ 
cal  effects.  Studies  on  the  cellular  level  show  increased  transport 
across  the  membrane  at  much  lower  exposure  times  or  number  of 
pulses.  Exposed  to  2000  pulses,  NG108  cells  exhibit  an  increase  in 
membrane  conductance,  but  only  allow  transmembrane  currents 
to  flow,  if  the  medium  is  positively  biased  with  respect  to  the  cell 
interior.  This  means  that  the  cell  membrane  behaves  like  a  rectify¬ 
ing  diode.  This  increase  in  membrane  conductance  is  a  nonthermal 
process,  since  the  temperature  rise  due  to  the  pulsing  is  negligible. 

Index  Terms — Cell  death,  current  rectifying,  permeabilization, 
pulsed  electric  field,  subnanosecond  pulses. 


I.  Introduction 

INTENSE  electrical  pulses  with  durations  of  microseconds 
to  milliseconds  are  known  to  cause  electroporation  [1]— [7] . 
The  pulses  charge  the  membrane  and  cause  an  increase  in  trans¬ 
membrane  voltage  from  250  mV  to  1  V  [8],  at  which  point 
membrane  permeabilization  or  pore  formation  sets  in.  This  pore 
formation  allows  delivery  of  drugs  or  genes  into  the  cell.  If  the 
pores  close  after  the  pulses,  the  process  is  known  as  reversible 
electroporation.  Direct  cell  killing  can  also  be  achieved  with 
pulses  of  the  same  duration  but  of  higher  amplitude,  referred  to 
as  irreversible  electroporation  [9] .  However,  when  shorter  pulses 
are  applied,  with  durations  in  the  submicrosecond  range,  inter¬ 
nal  effects  have  been  observed  without  permanent  damage  to 
the  cell  membrane  (intracellular  electromanipulation)  [10].  The 
effects  include  internal  calcium  bursts  [11],  [12],  permeabiliza¬ 
tion  of  internal  granules  [13],  and  control  of  cell  death  through 
apoptosis.  The  transition  from  electroporation  to  intracellular 
electromanipulation  occurs  for  pulse  durations,  which  are  on  the 
order  of  the  charging  time  constant  of  a  cell  membrane,  which 
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is  approximately  100  ns  for  mammalian  cells.  Yet,  the  duration 
of  the  pulses  is  not  the  only  criterion  for  the  onset  of  electro¬ 
poration  or  intracellular  electroporation.  The  peak  electric  field 
strength  is  important  too.  For  example,  nanosecond  pulses  with 
high  intensity  and  fast  rise  time  can  also  cause  membrane  effects, 
such  as  membrane  permeabilization  [14],  [15],  nanometer- sized 
pore  formation  [16],  and  phosphatidylserine  (PS)  extemaliza- 
tion  [17].  In  fact,  increasing  the  number  of  nanosecond  pulses 
can  cause  the  permeabilization  of  the  membranes  of  endocytotic 
vesicles  as  well  as  the  plasma  membrane  [18].  Nonetheless,  one 
could  still  expect  that  the  shorter  the  pulse  duration  is,  the  higher 
the  likelihood  that  the  pulses  can  interact  with  the  cell  inter¬ 
nal  organelles  is.  For  subnanosecond  pulses,  the  probability  of 
penetrating  into  the  interior  of  the  cell  is  even  higher  than  for 
nanosecond  pulses  [19].  In  this  case,  the  electric  field  distribu¬ 
tion  is  determined  by  the  dielectric  permittivity  rather  than  the 
resistivity  of  cell  components  [20]. 

There  are  two  reasons  to  enter  the  subnanosecond  range  of 
pulse  electric  fields.  The  first  is  that  for  extremely  short  pulses, 
the  electric  field  acts  directly  on  membrane  proteins,  rather  than 
causing  charging  of  the  membrane,  and,  if  sufficiently  strong, 
may  cause  direct  and  instant  conformational  changes.  A  model 
analysis  [21]  has  shown  that  subnanosecond  pulses  can  perforate 
a  cell  membrane.  In  addition,  subnanosecond  pulses  (800  ps) 
have  been  shown  to  cause  B16  cells  to  increase  trypan  blue 
uptake  [20],  which  is  a  clear  indication  of  the  disruption  of 
membrane  integrity.  The  second  reason  is  that  subnanosecond 
pulses  may  allow  us  to  use  an  impulse  antenna  [22],  [23]  to 
deliver  such  pulses  for  noninvasive  treatment  of  deep  lesions. 
The  preferred  pulse  duration  for  subnanosecond  pulses  is  in  the 
range  of  100-200  ps,  making  it  possible  to  focus  the  radiation 
on  the  target  efficiently  and  produce  a  focal  spot  of  1-cm  size  in 
the  tissue. 

In  this  paper,  we  present  a  study  of  cell  responses  to  200  ps 
(70%  peak  width).  The  biological  cell  response  to  such  short 
pulses  has  not  been  investigated  to  the  best  of  our  knowledge. 
Cell  survival  for  different  pulse  numbers  was  assessed  with  the 
water  soluble  tetrazolium  (WST)  assay,  18  h  after  exposure.  The 
membrane  transport  process  immediately  after  exposure  was 
studied  using  patch  clamp.  Our  goal  was  to  assess  the  primary 
response  of  the  cell  membrane  within  1  min  after  the  exposure, 
as  well  as  to  evaluate  the  cell  death  18  h  after  the  exposure. 

II.  Materials  and  Methods 

Two  types  of  cells  were  studied:  B16.F10  murine  melanoma 
cells  (ATCC  CRL-6322)  and  NG108  neuroblastoma-glioma  hy¬ 
brid  cells  (ATCC  HB-12317).  B16  cells  were  used  for  the  cell 
survival  study.  NG108  cells  have  voltage-gated  (VG)  sodium 
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Fig.  1 .  Subnanosecond  exposure  system,  (a)  Cuvette  with  a  loading  volume  of 
6  n L.  The  electrodes  are  70-/im-thick  stainless  steel  foils,  (b)  Cuvette  was  placed 
at  the  end  of  a  coax-cable  (RG  217).  A  conical  resistor  array  was  connected  in 
parallel  with  the  cuvette  for  matching  the  50-0  coaxial,  subnanosecond  pulse 
delivery  cable. 

channel  current  that  is  sensitive  to  the  electric  pulses  and,  there¬ 
fore,  were  used  for  studying  the  change  of  transmembrane  con¬ 
ductance  due  to  the  electric  pulses.  The  pulse  conditions  for 
B16  cell  viability  study  were  also  used  on  NG108  cells  to  as¬ 
sess  the  cell  death.  The  cells  were  maintained  as  monolayers  in 
90%  McCoy’s  medium  (B16)  or  Dulbecco’s  Modified  Eagle’s 
Medium  without  sodium  pyruvate  with  10%  fetal  bovine  serum 
(NG108).  The  media  ingredients  were  purchased  from  Gibco 
(Grand  Island,  NY).  B16  cells  were  removed  from  flasks  us¬ 
ing  a  cell  nonenzy mafic  dissociation  solution  (Sigma,  St.  Louis, 
MO)  [24]. 

To  evaluate  the  effects  that  subnanosecond  pulses  have  on  bi¬ 
ological  cells,  we  used  a  specially  designed  cuvette  with  plane- 
parallel  electrodes,  2  mm  apart,  for  cell  exposure.  The  geometry 
of  the  cuvette  is  shown  in  Fig.  1(a).  The  electrodes  are  stainless 
steel  foil  with  a  thickness  of  70  /im.  The  cuvette  was  loaded 
with  6-pL  cell  buffer.  The  reason  for  choosing  such  a  small 
volume  was  to  minimize  the  capacitance  of  the  cuvette  in  order 
to  ensure  that  the  fast  rise  time  of  the  pulse  will  not  be  slowed. 
The  total  capacitance  is  calculated  as  1 .4  pF,  which  allows  the 
shortest  rise  time  of  the  pulse  to  be  70  ps  for  a  50-U  coaxial 
delivery  system.  The  pulses  in  our  experiment  have  a  rise  time 
of  200  ps  (14%-92%  pulse  amplitude)  and,  therefore,  retained 
their  fidelity  across  the  cuvette.  The  second  reason  was  that  the 
cuvette  had  to  be  small  in  order  to  be  placed  in  a  coaxial  well 
[see  Fig.  1(b)].  The  cuvette  was  placed  in  parallel  with  a  conical 
resistor  array  that  was  connected  to  the  delivery  cable  from  the 
pulse  source.  The  conical  resistor  array  (total  value:  48.9  Cl) 
allows  for  a  continuous  potential  distribution  from  the  delivery 
cable  and  terminates  the  power  flow  without  creating  strong  re¬ 
flections.  For  a  200-ps  pulse,  the  He  bandwidth  is  approximately 
3  GHz.  We  used  an  HP  network  analyzer  (HP8753E)  to  measure 
the  returned  loss  for  frequencies  from  30  kHz  to  3  GHz  and  it 
was  found  to  be  between  —30  dB  and  —  lldB  in  this  range; 
therefore,  the  cuvette  was  well  matched  to  the  pulse  delivery 
cable,  and  the  reflection  was  minimized.  The  subnanosecond 
pulses,  with  an  amplitude  of  5  kV  were  generated  by  a  pulse 
generator  FPG  5-P  (FID  GmbH,  Germany);  therefore,  the  elec¬ 
tric  field  was  25  kV/cm.  This  electric  field  is  far  less  than  the 


breakdown  field  of  4  MV/cm  in  water  and  75  kV/cm  in  air 
for  20-ns-long  electric  pulses  [25],  [26].  For  subnanosecond 
pulses,  the  breakdown  field  is  even  higher  in  both  air  and  water 
than  for  nanosecond  pulses.  Therefore,  we  eliminate  the  pos¬ 
sibility  of  electric  breakdown  inside  the  cuvette.  During  puls¬ 
ing,  the  temperature  was  monitored  by  a  fiberoptic  temperature 
sensor  (T1  C- 11000  A,  Neoptix,  Canada)  with  a  tip  diameter 
^300  fi m.  The  temperature  probe  was  positioned  on  a  microm¬ 
eter  stage  and  was  dipped  into  the  cell  buffer  from  the  opening 
of  the  cuvette.  The  sampling  frequency  was  2  Hz. 

A  separate  electrode  configuration  was  used  to  study  the  trans¬ 
membrane  conductance  with  the  patch  clamp  technology.  The 
pulses  were  delivered  to  cells  on  a  cover  slip  on  a  microscope 
stage  by  a  mini-coax  cable  (300  /im  outer  diameter,  and  70 
/im  diameter  of  the  inner  conductor).  The  mini-coax  has  its 
end  opened  with  the  inner  conductor  exposed  (exposed  length: 
350  /im).  With  the  help  of  a  robotic  manipulator  (MP-225, 
Sutter),  the  electrode  was  positioned  above  the  cover  slip  sur¬ 
face,  and  the  angle  between  the  cover  slip  and  the  electrode 
was  20°.  The  leakage  field  from  the  inner  conductor  to  the 
outer  conductor  provides  the  exposure  environment  for  the  cells. 
The  electric  field  at  the  cell  location  on  the  cover  slip  between 
the  electrodes  was  calculated  by  3-D  simulations  with  a  finite- 
element  Laplace  equation  solver  Amaze  3-D  (Field  Precision, 
Albuquerque,  NM).  The  electric  field  was  calculated  to  be  20- 
40  kV/cm  in  the  exposure  area  of  200  /im  x  200  fi m  on  the 
cover  slip.  Membrane  currents  were  measured  before  exposure 
to  subnanosecond  pulses  and  at  10  and  20  s  after  it.  The  strong 
inhomogeneity  of  the  field  prevents  us  from  exposing  a  large 
number  of  cells  to  the  same  electric  field  intensity.  Therefore, 
this  pulse  system  was  not  used  for  cell  viability  study  using  the 
WST  assay. 

III.  Results  and  Discussions 
A.  Cell  Survival  After  Exposure  to  Subnanosecond  Pulses 

Our  previous  study  [20]  showed  that  for  the  800-ps  pulses 
and  an  electric  field  of  150  kV/cm,  approximately  18  000  pulses 
were  required  in  order  to  cause  50%  B-16  cell  death.  At  the 
present  state  of  antenna  technology,  such  electric  fields  cannot 
be  generated  in  tissues  by  means  of  wideband  antennas.  We, 
therefore,  set  the  electric  field  to  25  kV/cm.  Compared  to  the 
earlier  study  [20],  both  the  electric  field  and  the  pulse  dura¬ 
tion  were  decreased  considerably.  It  raised  the  question  of  how 
effective  these  short  pulses  can  be  in  eliminating  unwanted  cells. 

An  initial  assessment  of  B16  cell  death  for  the  200-ps  pulses 
at  25  kV/cm  was  performed  under  two  pulse  conditions,  1000 
pulses  (300  Hz)  and  10000  pulses  (300  Hz).  Using  a  WST 
assay  1 8  h  after  the  pulsing,  we  observed  no  significant  effect  of 
exposure.  The  results  were  comparable  to  the  control  (in  which 
cells  were  cultured  in  the  incubator  all  along)  and  the  sham 
exposure  (when  cells  were  loaded  into  the  exposure  cuvette,  but 
no  pulse  was  applied). 

In  order  to  achieve  significant  cell  death,  the  number  of  pulses 
has  to  be  increased  considerably.  The  repetition  rate  also  needs 
to  be  increased,  as  it  was  shown  [27]  that  increasing  the  pulse 
repetition  rate  caused  an  increase  of  cell  killing,  even  though  the 
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Fig.  2.  Temperature  increases  due  to  the  application  of  1.2  M,  1.8  M,  and 
2.4  M  pulses.  Here,  “M”  stands  for  a  million.  The  pulses  were  applied  in  the 
burst-mode  with  20, 30,  and  40  pulses/burst.  The  repetition  rate  in  each  burst  was 
1 0  kHz  and  each  burst  repeated  every  1 0  ms  ( 1 00  Hz) .  The  exposure  time  for  each 
of  these  three  exposure  conditions  was  10  min.  The  starting  temperature  was 
24°C±  1  °C. 


applied  voltage  was  lower  than  the  membrane  breakdown  volt¬ 
age.  We  have  increased  the  repetition  rate  to  10  kHz,  the  max¬ 
imum  rate  the  pulse  generator  allowed.  Thermal  effects,  which 
generally  can  be  neglected  for  nanosecond  and  subnanosecond 
pulses  at  low  repetition  rates  must  be  taken  into  account  at  high 
repetition  rates.  We  have  found  that  the  pulses  delivered  to  a 
cuvette  at  high  repetition  rates  caused  a  significant  increase  in 
the  temperature  of  the  cell  buffer.  For  a  10-min  exposure  at 
3  kHz,  a  total  of  1.8  million  pulses,  at  maximum  temperature 
of  37  °C,  was  measured.  Results  obtained  at  10  kHz,  showed 
a  much  higher  absolute  temperature  (>60  °C).  The  exposure 
at  the  increased  temperature  also  caused  evaporation  of  the  cell 
buffer  and  the  loss  of  volume  of  the  buffer,  seen  from  the  volume 
retrieved  through  a  micropipette  with  a  gel-loading  tip  (Fisher 
Scientific,  PA). 

Instead  of  applying  pulses  continuously,  we  used  a  burst  mode 
to  deliver  the  pulses.  Each  burst  consisted  of  20,  30,  or  40  pulses 
and  the  burst  repeated  itself  every  10  ms  (100  Hz).  The  pulses  in 
each  burst  had  a  repetition  frequency  of  10  kHz  (100  /is  between 
pulses).  For  20,  30,  and  40  pulses/burst,  the  total  pulse  numbers 
are  1.2,  1.8,  and  2.4  million  for  10-min  exposure.  The  tempera¬ 
ture  developments  are  shown  in  Fig.  2.  The  starting  temperature 
was  24  ±  1  °C.  Within  approximately  1  min,  the  temperature 
rose  and  became  stable.  The  corresponding  maximum  tempera¬ 
ture  was  measured  as  37  °C,  41  °C,  and  47  °C.  We  note  that  the 
temperature  in  Fig.  2  was  measured  when  the  temperature  probe 
tip  was  just  submerged  in  the  cell  buffer,  less  than  1  mm  deep. 
For  other  positions,  the  maximum  temperatures  decreased  to  34 
db  1  °C  for  20  pulses/burst,  37  =b  1  °Cfor30  pulses/burst,  and 
44  db  1  °C  for  40  pulses/burst,  respectively.  The  higher  tempera¬ 
ture  near  the  cell  buffer  surface  can  be  explained  by  the  electric 
field  enhancement  of  the  sharp  edge  electrodes  as  they  are  thin 
stainless  steel  (70-/im- thick)  foil.  The  field  near  the  surface  can 
be  estimated  at  10%  higher  than  the  average  field  in  the  gap.  On 
the  other  hand,  for  the  exposure  of  2.4  million  pulses,  the  energy 


Fig.  3.  Survival  of  B-16  cells  after  exposure  to  burst-mode  pulsing.  Pulse 
number,  0,  refers  to  sham  exposure.  Pulse  number  1.2  M,  1.8  M,  and  2.4  M 
corresponds  to  the  burst-mode  exposure  with  20,  30,  and  40  pulses/burst.  Here, 
“M”  stands  for  a  million.  The  error  bars  represent  the  standard  deviation  of  the 
mean  (n  =  10  for  each  pulsing  condition).  A  single  factor  ANOVA  analysis 
yielded  p  —  0.22  (1.2  M  pulses)  and  p  <  0.003  for  both  1.8  M  and  2  M  pulses. 


W  required  to  increase  the  temperature  of  the  6-fiL  cell  buffer 
for  a  temperature  increase,  AT  =  23  °C,  can  be  estimated  as 
0.6  J  (AW  =  CvpATV ,  Cv  =  4.2  J/g°C,  V=6  fi L).  The  overall 
energy  W  deposited  in  the  buffer  is  24  J,  calculated  from  W  = 
<jE2tNV,  where  a  =  1 .3  S/m,  E  =  25  kV/cm,  r  =  200  ps,  and 
N  =  2.4  million.  It  is  obvious  that  most  of  the  energy  was  dis¬ 
sipated  through  the  electrical  contact.  Fig.  3  shows  the  viability 
of  B16  cells  18  h  after  exposure.  For  exposure  to  1.2  million 
pulses  at  10  kHz  and  20  pulses/burst,  Fig.  3  indicates  a  slight  in¬ 
crease  of  cell  survival.  This,  however,  should  not  be  interpreted, 
as  a  higher  number  of  cells  proliferated  compared  to  the  sham. 
Instead,  this  tells  us  the  cells  increased  their  mitochondrial  ac¬ 
tivity  and  metabolic  rate  due  to  pulsing.  As  the  pulse  number 
increased,  the  number  of  cells  that  survived  pulsing  decreased. 
Approximately  20%  cell  death  was  observed  for  1.8  million 
pulses  (10  kHz  and  30  pulses/burst),  and  90%  cell  death  for  2.4 
million  pulses  (10  kHz  and  40  pulses/burst).  It  should  be  noted 
that  WST  is  a  colorimetric  assay,  which  measures  the  activity  of 
mitochondrial  enzymes.  It  is  likely  that  subnanosecond  pulses 
penetrate  into  the  cell  interior  and  affect  the  metabolic  activity 
directly,  and  it  is  possible  that  the  effect  could  occur  even  in 
viable  cells.  However,  we  observed  that  most  of  the  cells  treated 
with  2.4  million  pulses  were  isolated  and  detached  from  the 
culture  dish.  In  addition,  when  cell  viability  was  assessed  by 
adding  trypan  blue,  there  was  a  clear  indication  that  most  of  the 
cells  were  dead.  In  contrast,  live  cells  in  sham  conditions,  which 
served  as  a  control,  were  confluent  and  remained  attached  to  the 
culture  dish  (as  well  as  appearing  viable  by  WST).  Therefore, 
we  think  the  result  obtained  by  the  WST  assay  is  a  valid  means 
for  determining  cell  death.  In  order  to  better  understand  whether 
the  temperature  rise  was  the  cause  of  the  observed  cell  death, 
we  incubated  cells  in  a  water  bath  at  various  temperatures  for  10 
min.  Fig.  4(a)  shows  the  temperature  profile  during  the  10-min 
water  bath  treatment  at  42  °C  when  the  cells  were  loaded  into 
a  test  tube.  The  temperature  in  the  tube  rose  to  the  water  bath 
temperature  in  approximate  50  s,  which  is  a  similar  trend  as  in 
the  pulsed  conditions.  The  survival  results  were  again  obtained 
from  the  WST  assay  18  h  after  exposure  and  shown  in  Fig.  4(b). 
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Fig.  4.  Cell  viability  after  treatment  in  water  baths  at  different  temperatures, 
(a)  Temperature  profile  for  10-min  treatment  was  measured  in  a  plastic  test  tube 
loaded  with  15-/iL  cell  buffer,  (b)  B-16  cell  survival  after  exposure  to  heating  in 
water  baths  at  various  temperatures  up  to  15  min.  The  error  bars  represent  the 
standard  deviation  of  the  mean  (n  =  4).  The  time  does  not  exactly  correspond 
to  the  treatment  time  at  the  water  bath  temperature,  as  the  temperature  in  the 
test  tube  took  approximately  50  s  to  rise.  At  45  °C,  the  15-min  data  are  not 
significantly  different  from  the  10-min  data  with  p  =  0.0892. 

The  cells  that  were  exposed  to  42  °C  were  alive,  even  after  the 
10-min  exposure.  However,  when  the  water  bath  temperature 
was  raised  to  47  °C,  almost  all  the  cells  did  not  survive.  Recall 
that  the  temperature  in  the  cuvette  under  the  pulse  condition  of 
40  pulses  per  burst  was  44  °C  -47  °C.  For  cells  treated  in  the 
water  bath  at  44  °C  -47  °C,  we  observe  that  cell  viability  could 
vary  in  the  range  of  20%-80%  [see  Fig.  4(b)].  Evidently,  the 
heating  contributed  to  the  cell  death.  We,  therefore,  conclude 
that  the  cell  death  under  the  pulse  condition  of  40  pulses  per 
burst  in  Fig.  3  was  caused  by  a  combination  of  heating  and  the 
electric  field  effect.  However,  for  the  pulse  condition  10  kHz  and 
30  pulses/burst,  the  heating  effect  was  minimal,  as  the  absolute 
temperature  was  37  °C-41  °C,  which  was  insufficient  to  cause 
cell  death,  as  shown  in  Fig.  4(b). 

B.  Membrane  Permeabilization  Caused  by  Subnanosecond 
Pulses 

We  have  observed  that  subnanosecond  pulses  caused  cell 
death  with  a  significant  temperature  rise  due  to  a  high-repetition 
rate  and  high  pulse  numbers.  Exposure  conditions  with  a  re¬ 
duced  pulse  number  and  a  negligible  temperature  rise  appear 
not  to  induce  any  noticeable  cell  death.  There  are,  however, 
other  cell  responses,  such  as  a  change  in  cell  membrane  con¬ 


ductance,  which  may  occur  due  to  the  applied  electric  field, 
especially  at  the  time  of  the  exposure  or  right  after  the  exposure. 
Under  normal  physiological  conditions,  the  transmembrane  con¬ 
ductance  is  primarily  determined  by  passive  transport  of  ions 
through  channels.  Ion  channels  are  pore-forming  proteins  and 
they  help  ions  move  across  the  plasma  membrane  according  to 
electrochemical  gradients.  The  most  important  are  two  types: 
ligand-gated  ion  channels  and  VG  channels. 

In  our  study,  we  employed  whole-cell  patch  clamp  to  explore 
the  effect  of  subnanosecond  electric  pulses  on  the  VG  sodium 
current  (VG  ?Na)  and  the  leak  current  (/Leak)  in  neuroblastoma 
cells  (NG-108).  All  experiments  were  carried  out  at  room  tem¬ 
perature  (^24  °C).  Whole-cell  patch-clamp  configuration  was 
established  2-3  min  prior  to  200-ps  electric  pulse  (psEP)  ap¬ 
plication,  to  allow  enough  time  for  dialysis  of  the  cytoplasm 
with  the  pipette  solution.  Whole-cell  currents  were  probed  by 
stepping  the  membrane  potential  from  —100  mV  to  +30  mV, 
in  10-mV  steps  from  the  holding  potential  of  —80  mV.  The  cur¬ 
rents  were  recorded  10  s  before  psEPs,  and  10  s  and  30  s  after 
the  exposure.  Shown  in  Fig.  5  are  the  current-voltage  curves  of 
VG  /ya  and  /Leak  in  cells  subjected  to  0  (sham,  n  =  4),  1000 
pulses  (n  =  6),  and  2000  pulses  (n  =  9)  applied  at  the  rate  of 
0.5  kHz.  VG  /Na  began  to  activate  at  a  potential  of  —40  mV, 
reached  peak  amplitude  near  10  mV,  and  then  declined  as  the 
voltage  approached  its  equilibrium  potential.  No  significant 
change  of  the  current-voltage  characteristics  of  VG  I^a  was 
observed  after  1000  pulses  and  even  2000  pulses  [see  Fig.  5(b) 
and  (c)].  The  /Leak  in  NG-108  cells  after  the  exposure  of  1000 
psEP  also  did  not  change  [see  Fig.  5(e)].  However,  /Leak  in  2000 
psEP-treated  cells  demonstrated  a  slight  increase  in  voltage- 
sensitive  inward  current  [see  Fig.  5(f)].  It  can  be  seen  that  10  s 
after  exposure,  the  /Leak  in  the  negative  voltage  had  increased 
11.1  ±  3pA. 

This  increase  was  still  observable  30  s  after  the  exposure,  but 
a  trend  of  recovery,  i.e.,  the  current  level  returning  to  the  control 
level,  could  be  seen.  In  the  meantime,  the  current  flowing  out 
of  the  cell  did  not  change.  This  unidirectional  current  was  most 
likely  due  to  the  influx  of  sodium  ions,  the  major  cation  pre¬ 
sented  in  the  cell  bath.  Such  unidirectional  current  flow  has  also 
been  observed  in  the  cell  membrane  after  exposure  to  nanosec¬ 
ond  electric  pulses  [14]  and  was  found  to  be  caused  by  pores  of 
nanometer  size  formed  in  the  cell  membrane.  When  nanosecond 
pulses  are  used  as  stimuli,  this  process  is  referred  to  as  nanopora- 
tion.  Interestingly,  the  current-rectifying  characteristic  of  these 
pores  is  similar  to  that  of  synthetic  nanopores  in  polymer  foils, 
reported  to  be  an  asymmetrical  conical  shape  [28]. 

These  pores,  as  they  were  created  by  a  relatively  small  num¬ 
ber  of  pulses,  were  unlikely  to  cause  cell  death,  since  the  cell 
viability  assessed  by  WST  assay  18  h  after  the  exposure  of 
2.4  million  pulses  was  found  to  be  34%  ±  2%  (temperature: 
42  °C-45  °C  during  exposure).  Cell  viability  after  exposure  to 
1.8  million  pulses  was  70%  (see  Fig.  6).  Once  the  pores  are 
formed,  can  they  develop  into  bigger  ones,  if  more  pulses  are 
applied?  Can  these  pores  become  large  enough  to  allow  trypan 
blue  molecules  to  enter  the  cell? 

One  hour  after  exposure  (1. 2-2.4  million  pulses),  the  cells 
did  not  show  more  trypan  blue  uptake  than  the  control.  This 
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Fig.  6.  Survival  of  NG108  cells  after  exposure  to  burst-mode  pulsing  using 
the  exposure  system  shown  in  Fig.  1.  Pulse  number,  0,  refers  to  sham  exposure. 
Pulse  number  1.2  M,  1.8  M,  and  2.4  M  corresponds  to  the  burst-mode  exposure 
with  20,  30,  and  40  pulses/burst.  Here,  “M”  stands  for  a  million.  The  error 
bars  represent  the  standard  deviation  of  the  mean  (n  =  4  for  each  pulsing 
condition).  The  maximum  temperatures  were  34  °C  (20  pulses/burst),  38  °C 
(30  pulses/burst),  and  45  °C  (40  pulses/burst). 
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Fig.  5.  Current-voltage  characteristics  of  NG108  cell  membrane  before  and 
after  exposure  to  subnanosecond  pulses.  Subnanoseocnd  electric  pulses  (psEPs) 
had  no  effect  on  VG  7n&  ( left  column )  but  increased  /Leak  (right  column ).  The 
currents  were  recorded  10  s  before  psEPs  (control),  and  10  s  and  30  s  after 
exposure.  Peak  amplitude  of  the  fast  inactivating  VG  [see  (a)-(c)]  was 
measured  as  a  negative  peak  during  the  first  15  ms  of  the  voltage  steps.  /Leak 
[see  (d)-(f)]  was  measured  at  a  mean  current  value  20-70  ms  after  stepping, 
when  VG  /Na  completely  inactivated.  Data  are  shown  as  the  mean  ±  SE  with  n 
=  4  for  sham,  n  =  6  for 1000  pulses,  and  n  =  9  for  2000  pulses. 


indicates  that  even  though  there  might  be  pores  formed  in  the 
cell  membrane,  their  size  is  not  large  enough  to  allow  trypan  blue 
to  enter  the  cell.  Another  possibility  is  that  these  pores,  even  if 
they  are  formed  and  become  larger  in  size,  might  reseal  in  less 
than  1  h,  as  has  been  observed  with  800-ps-long  pulses  [20]. 
One  other  possibility  is  the  pores  may  open  to  their  largest 
size  during  pulsing,  which  may  allow  trypan  blue  molecules  to 
enter  the  cell  during  the  exposure.  To  test  this  ideas,  we  added 
trypan  blue  into  the  cell  buffer  before  pulsing,  and  then  applied 
the  pulses  to  the  cells  in  the  buffer  under  the  same  conditions. 
Still,  we  were  not  able  to  observe  any  more  trypan  blue  uptake 
compared  to  the  sham. 

The  absence  of  large  pores  in  the  cell  membrane  was  also 
shown  in  B16  cells.  We  added  bleomycin  (Sigma,  molecu¬ 
lar  weight  1415.6Da),  a  chemotherapeutic  agent,  to  the  cell 


Fig.  7.  B-16  cell  survival  for  various  drug  (bleomycin)  concentrations.  The 

cell  survival  for  bleomycin  only  and  bleomycin  delivery  with  microsecond 
pulses  (100  /is)  serve  as  sham  and  positive  control.  The  number  of  microsecond 
pusle  applied  was  8.  The  number  of  subnanosecond  pulses  was  0.5  million.  The 
error  bars  represent  the  standard  deviation  of  the  mean  (n  =  3). 

buffer  [29].  We  then  applied  the  pulses  to  the  cell  buffer  to  see 
if  bleomycin  could  enter  the  cell.  The  pulse  number  in  this  case 
was  500  000,  still  not  sufficient  to  cause  cell  death.  The  cells 
were  again  cultured  for  18  h  prior  to  the  WST  assay.  The  sur¬ 
vival  result  is  shown  in  Fig.  7.  It  can  be  seen  that  the  survival  of 
the  cells  in  the  buffer  with  drug  alone  was  reduced  as  the  dose 
of  bleomycin  increased.  This  is  an  expected  result.  For  cells 
treated  with  subnanosecond  pulses,  however,  the  survival  was 
approximately  the  same  as  that  of  adding  bleomycin  only.  Ob¬ 
viously,  the  added  drug  did  not  cause  an  increase  in  cell  death. 
In  contrast,  for  electroporation  pulses  (100  fi s),  which  served 
as  a  positive  control,  the  cell  survival  was  decreased  at  a  much 
lower  bleomycin  concentration  (2  x  10-8  mol)  in  comparison 
to  the  toxic  threshold  (2  x  10-5  mol)  of  the  drug  alone.  The 
electroporation  pulses  increased  the  rate  of  cell  death  due  to 
uptake  of  bleomycin  compared  to  using  the  drug  alone  or  the 
use  of  subnanosecond  pulses. 

Both  the  trypan  blue  staining  study  and  bleomycin  study  sug¬ 
gest  that  the  nanopores,  once  they  are  formed,  may  not  become 
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larger  in  size  even  though  more  pulses  are  applied.  This  may  be 
due  to  the  fact  that  the  mechanism  of  such  short  pulses  inter¬ 
acting  with  cell  membrane  is  different  than  that  of  long  pulses. 
For  subnanosecond  pulses,  displacement  current  dominates  the 
total  current  and  the  conduction  current  flowing  through  the 
pores  is  much  less.  Even  though  lower  resistance  nanopores 
formed,  the  majority  of  the  current  still  passes  through  the  rest 
of  the  cell  membrane  as  displacement  current.  Conversely,  for 
long  electroporation  pulses,  the  pore  size  might  grow  larger,  if 
more  pulses  are  applied.  This  expansion  can  be  qualitatively  un¬ 
derstood  by  a  positive-feedback  mechanism.  As  the  conduction 
current  flows  through  the  cell  membrane,  the  current  mostly 
chooses  the  high-conductivity  path  via  pores.  This  results  in 
high  current  density  in  these  pores.  The  energy  density  within 
these  narrow  channels  is,  therefore,  significantly  higher  than 
the  nearby  intact  membrane.  The  temperature  in  the  current 
path  can,  therefore,  be  further  enhanced,  and  results  in  higher 
conductivity.  Such  positive-feedback  mechanism  may  allow  the 
pores  to  increase  in  diameter.  To  summarize,  the  pore  formed 
in  the  subnanosecond  range  is  most  likely  caused  by  the  direct 
electric-field-induced  conformational  change,  a  strictly  nonther- 
mal  process. 

IV.  Conclusion 

Subnanosecond  pulses,  much  shorter  than  the  cell  membrane 
charging  time  constant,  were  found  to  cause  cell  death  and  a 
change  in  membrane  permeability.  For  the  electric  pulse  con¬ 
dition,  200  ps,  25  kV/cm,  using  1.8  million  pulses,  cell  death 
(^20%)  was  observed  18  h  after  exposure.  The  cell  death  was 
caused  by  temperature-assisted  electric  field  effects.  A  higher 
number  of  pulses  (2.4  million)  resulted  in  a  greater  temperature 
rise,  and  90%  cell  death  was  observed.  In  either  case,  the  cell 
membrane  was  found  intact  1  h  after  exposure  when  tested  by 
trypan  blue. 

While  the  pulses  with  durations  shorter  than  the  charging  time 
of  a  cell  membrane  (^100  ns  for  mammalian  cells)  have  a  higher 
probability  of  acting  upon  intracellular  structures,  this  proposi¬ 
tion  should  not  create  the  impression  that  cell  membranes  will 
be  completely  exempt  from  the  effects  of  the  subnanosecond 
pulses.  Through  the  patch  clamp  study  on  NG108  cells,  we  have 
found  that  the  cell  membrane  increased  its  conductance  in  re¬ 
sponse  to  subnanosecond  electric  pulses  (electric  field  between 
20-40  kV/cm,  200  ps).  Also,  the  increase  in  the  conductance 
lasted  for  at  least  30  s  after  pulse  application.  This  increased 
conductance  occurred  at  much  lower  number  of  pulses  than  re¬ 
quired  for  causing  significant  cell  death,  which  suggests  that  this 
process  is  nonthermal.  In  addition,  this  process  is  not  necessar¬ 
ily  associated  with  cell  death.  More  interestingly,  the  increased 
conductance  only  allowed  a  unidirectional  current  flow  from  the 
extracellular  to  the  intracellular  medium.  This  unidirectionality 
of  current  flow  exhibited  the  same  features  as  nanopores  formed 
on  the  cell  membrane  after  being  exposed  to  nanosecond  pulses. 
We  note  that  if  nanopores  form  in  the  cell  membrane,  the  pores 
may  not  grow  larger  in  size,  even  though  more  pulses  are  applied. 
These  pores  are  not  large  enough  to  allow  either  trypan  blue  or 
bleomycin  molecules  to  enter  the  cell  interior.  However,  this 


scenario  may  change,  if  one  applies  higher  electric  fields,  which 
may  be  sufficient  to  cause  stronger  conformational  changes  in 
the  cell  membrane. 

In  general,  when  using  subnanosecond  electric  pulses  for 
therapeutic  applications  in  eliminating  unwanted  cells,  we  con¬ 
clude  that  an  electric  field  of  at  least  20  kV/cm  is  required. 
The  cell  death  is  due  to  both  electric-field  and  thermal  effects. 
Delivering  pulses  of  higher  electric  field  intensity  at  a  lower 
repetition  rate  may  allow  us  to  separate  the  effects,  which  jus¬ 
tifies  further  study.  On  the  other  hand,  subnanosecond  pulses 
offer  a  potential  treatment  modality  that  can  modify  cell  mem¬ 
brane  properties  or  alter  the  physiological  conditions  of  cells  to 
induce  a  desired  response  instead  of  causing  cell  death.  In  this 
case,  much  lower  numbers  of  pulses  are  needed  to  obtain  this 
effect.  Such  applications  can  include  enhancement  of  wound 
healing,  or  hemostasis  (blood  clotting),  by  stimulating  calcium 
influx  into  platelets  [30] . 
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Abstract — Nanosecond  electrical  pulses  have  been  successfully 
used  to  treat  melanoma  tumors  by  using  needle  arrays  as  pulse 
delivery  systems.  Reducing  the  pulse  duration  of  intense  electric 
field  pulses  from  nanoseconds  into  the  subnanosecond  range  will 
allow  us  to  use  wideband  antennas  to  deliver  the  electromagnetic 
fields  into  tissue  with  a  spatial  resolution  in  the  centimeter  range. 
To  explore  the  biological  effect  of  intense  subnanosecond  pulses, 
we  have  developed  a  generator  that  provides  voltage  pulses  of 
160  kV  amplitude,  200  ps  rise  time,  and  800  ps  pulse  width.  The 
pulses  are  delivered  to  a  cylindrical  Teflon  chamber  with  polished 
flat  electrodes  at  either  end.  The  distance  between  the  electrodes 
is  variable  and  allows  us  to  generate  electric  fields  of  up  to 
1  MV/cm  in  cell  suspensions.  The  pulses  have  been  applied  to  B16 
(murine  melanoma)  cells,  and  the  plasma  membrane  integrity  was 
studied  by  means  of  trypan  blue  exclusion.  For  pulse  amplitudes 
of  550  kV/cm,  approximately  50%  of  the  cells  took  up  trypan 
blue  right  after  pulsing,  whereas  only  20%  were  taking  it  up  after 
1  h.  This  indicates  that  the  plasma  membrane  in  a  majority  of 
the  cells  affected  by  the  pulses  recovers  with  a  time  constant  of 
about  1  h.  The  cells  that  show  trypan  blue  uptake  after  this  time 
suffer  cell  death  through  apoptosis.  Evaluation  of  the  experimental 
results  and  molecular  dynamics  modeling  results  indicate  that 
with  a  pulse  duration  of  800  ps,  membrane  charging  and  nanopore 
formation  are  the  dominant  bioelectric  effects  on  B16  cells.  This 
information  has  been  used  in  a  continuum  model  to  estimate 
the  increase  in  membrane  permeability  and,  consequently,  the 
increase  in  pore  size  caused  by  repetitive  pulsing. 

Index  Terms — Apoptosis,  biological  cells,  cell  membranes, 
nanopores,  picosecond  pulses,  pulse  power,  wideband  antenna. 

I.  Introduction 

THE  EFFECT  of  high-intensity  nanosecond  electrical 
pulses  on  biological  cells  and  tissue  has  received  con¬ 
siderable  attention  by  the  scientific  community,  particularly 
in  the  past  five  years  and  has  led  to  the  establishment  of  a 
new  research  field:  bioelectrics  [1],  [2].  In  this  time  domain, 
the  rise  time  of  the  pulse  is  faster  than  the  charging  time  of  the 
plasma  membrane  of  most  mammalian  cells,  meaning  that  the 
field  will  pass  through  the  membrane  into  the  cytoplasm.  This 
effect  can  be  understood  qualitatively  by  considering  the  cell  as 
an  electrical  circuit,  describing  the  various  cell  membranes  by 
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Fig.  1.  (a)  Structure  of  a  biological  cell  (as  would  be  seen  with  a  light 

microscope),  (b)  Double- shell  model  of  a  biological  cell  and  superimposed 
equivalent  circuit  of  the  cell  [2]. 


their  capacitances  and  the  cytoplasm,  which  they  enclose,  by 
its  resistance.  Fig.  1(a)  shows  the  cross  section  of  a  mammalian 
cell,  with  the  only  membrane-bound  substructure  shown  being 
the  nucleus.  The  cytoplasm,  which  fills  much  of  the  cell,  con¬ 
tains  dissolved  proteins,  electrolytes,  and  glucose  and  is  mod¬ 
erately  conductive,  as  are  the  nucleoplasm  and  the  cytoplasm 
in  organelles.  On  the  other  hand,  the  membranes  that  surround 
the  cell  and  subcellular  structures  have  low  conductivity.  We 
can  therefore  think  of  the  cell  as  a  conductor  surrounded  by 
an  ideally  insulating  envelope,  containing  substructures  with 
similar  properties.  The  equivalent  circuit  of  such  a  cell  (which 
is  considered  as  spherical  for  modeling  purposes)  with  one 
substructure,  the  nucleus,  is  shown  in  Fig.  1(b). 

If  dc  electric  fields  or  pulses  of  long  duration  (compared 
to  the  charging  time  of  the  capacitor  formed  by  the  outer 
membrane)  are  applied,  eventually,  only  the  outer  membrane 
will  be  charged;  the  electric  field  generated  across  subcellular 
membranes  during  charging  will  be  zero  for  an  ideal,  fully 
insulating  outer  membrane.  However,  during  the  charging  time 
of  the  outer  membrane,  we  will  also  expect  potential  differences 
to  be  generated  across  subcellular  membranes,  an  effect  that 
will  be  stronger  if  the  pulse  rise  time  is  shorter.  Such  charging 
times  are  in  the  submicrosecond  range  for  human  cells. 

If  the  field  is  sufficiently  large,  it  can  have  pronounced 
effects  on  intracellular  organelles.  High-voltage  pulses  that 
are  nanoseconds  to  hundreds  of  nanoseconds  long  have  been 
shown  to  penetrate  into  living  cells  to  permeabilize  intracellular 
organelles  [3],  [4]  and  release  Ca2+  from  the  endoplasmic 
reticulum  [5]— [7] .  They  provide  a  new  approach  for  physically 
targeting  intracellular  organelles  with  many  applications  in¬ 
cluding  activation  of  platelets  and  release  of  growth  factors  for 
accelerated  wound  healing  [8]  and  precise  control  of  apoptosis 
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[9]— [12] .  Most  recently  it  has  been  shown  that  such  pulsed 
electric  fields  cause  shrinkage  and  even  complete  elimination 
of  melanoma  tumors  [13]. 

Whereas  the  emphasis  on  bioelectric  applications  has  been 
on  nanosecond  pulse  effects,  there  are  two  good  reasons  to  ex¬ 
plore  the  effect  of  subnanosecond  pulses  on  biological  cells  and 
tissue.  The  first  is  that  for  extremely  short  pulses,  the  dielectric 
properties,  rather  than  the  resistive  characteristics  of  the  media, 
determine  the  electric  field  distribution.  In  the  equivalent  circuit 
shown  in  Fig.  1(b),  the  conductance  of  the  membranes  is  as¬ 
sumed  to  be  zero,  and  the  capacitive  components  of  cytoplasm 
and  nucleoplasm  are  neglected.  These  are  assumptions  that  are 
only  valid  for  a  pulse  duration  that  is  long  compared  to  the 
dielectric  relaxation  time  of  the  cytoplasm  {ecv/crcv),  where 
£cp  is  the  permittivity  of  the  cytoplasm  and  <rcp  its  conductivity. 
Only  then  can  the  capacitive  term  in  the  cytoplasm  impedance 
be  neglected  compared  to  the  resistive  term.  Based  on  the  data 
listed  in  [14],  this  does  not  hold  anymore  for  white  blood  cells 
when  they  are  exposed  to  pulse  durations  short  compared  to  one 
nanosecond.  The  electric  fields  in  the  various  parts  of  the  cell 
are  then  defined  by  their  permittivity,  rather  than  conductivity. 
For  a  membrane  with  a  relative  dielectric  constant  of  8,  the 
electric  field  in  the  membrane  is  ten  times  higher  than  the 
electric  field  in  the  adjacent  cytoplasm,  which  has  a  dielectric 
constant  of  80.  The  electric  field  then  acts  directly  on  membrane 
proteins,  rather  than  causing  charging  of  the  membrane,  and,  if 
sufficiently  strong,  can  cause  direct  and  instant  conformational 
changes,  such  as  voltage  gating. 

Second,  besides  providing  the  opportunity  to  enter  a  new 
domain  of  electric  field-cell  interactions,  subnanosecond  pulses 
will  ultimately  allow  medical  applications  for  delivery  of  pulsed 
electric  fields  without  invasive  electrodes,  using  antennas  in¬ 
stead.  A  possible  configuration  that  allows  us  to  generate  very 
high  electric  fields  uses  a  focusing  antenna  [15]  and  may  be 
able  to  induce  apoptosis  in  tissue  without  using  needles  as  the 
delivery  system  for  electric  pulses  [13]. 

To  explore  effects  of  electrical  pulses  with  durations  right 
at  the  theoretical  threshold  where  membrane-charging  effects 
might  become  less  important,  and  dielectric  effects  could  al¬ 
ready  be  present,  we  have  applied  800-ps-long  pulses  with 
electric  field  amplitudes  of  up  to  1  MV/cm  to  melanoma  (B16) 
cells.  We  have  measured  the  membrane  integrity  of  these  cells 
after  pulsing  and  determined  the  fraction  of  dead  cells  and  the 
type  of  cell  death,  depending  on  electric  field  intensity  and 
number  of  pulses. 

II.  Experimental  Setup  and  Methods 

A  pulsed  power  system,  which  is  able  to  provide  subnanosec¬ 
ond  pulses  to  a  biological  load,  was  designed  and  built  [16], 
[17]  (Fig.  2).  It  is  driven  by  a  waveguide  Marx-Bank  [18]  that 
charges  a  10-cm-long  coaxial  cable  (RG  213).  The  cable  is 
discharged  through  a  high-pressure  spark  gap.  The  system  is 
able  to  generate  800-ps-long  voltage  pulses  with  an  amplitude 
of  160  kV  into  a  50-9  load.  Voltages  measured  at  the  exposure 
cell,  which  acts  as  a  high-impedance  load,  reach  values  of 
approximately  300  kV.  A  typical  voltage  pulse  shape  is  shown 
in  Fig.  3.  With  a  gap  distance  of  3  mm,  we  have  reached  electric 
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Fig.  2.  Block  diagram  of  ultrashort  pulse  generator  with  cross  section  of 
exposure  system  (upper  right)  [17]. 
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Fig.  3.  Voltage  measured  at  the  exposure  cuvette  [17]. 

fields  of  up  to  1  MV/cm.  In  spite  of  this  extremely  large  field,  no 
electrical  breakdown  was  observed.  This  is  in  line  with  results 
obtained  with  200  ns  pulses,  where  the  breakdown  field  for 
water  exceeded  values  of  1  MV/cm  in  a  pin-plate  electrode 
configuration  [19].  With  pulse  durations  reduced  by  two  orders 
of  magnitude  compared  to  those  used  in  the  water  breakdown 
experiments,  it  is  not  surprising  that  even  multimegavolts  per 
centimeter  fields  in  the  subnanosecond  range  would  not  lead 
to  electrical  breakdown.  However,  streamer  formation  starting 
from  microprotrusions  at  the  electrodes  cannot  be  excluded  and 
needs  to  be  studied  [20] . 

A  murine  melanoma  cell  line,  B16,  was  chosen  as  a  mam¬ 
malian  cell  model  to  determine  the  effects  of  these  subnanosec¬ 
ond  pulses  on  membrane  permeability  and  cell  death.  The  cells 
in  suspension  were  placed  between  two  plane  parallel  stainless 
steel  electrodes  (cuvette)  with  a  variable  gap  distance  d  (3.15, 
3.91,  5.35,  and  10.67  mm,  respectively).  The  cylindrical  (flat- 
end)  electrodes  had  a  diameter  of  2.75  mm.  Cells  loaded  into 
the  cuvette  were  exposed  to  multiple  800  ps  pulses.  The  electric 
field  in  the  suspension,  V/d,  with  V  being  the  pulse  voltage, 
was  varied  between  150  kV/cm  and  1  MV/cm.  Cells  loaded 
into  cuvettes  and  merely  placed  into  the  circuit  without  being 
pulsed  were  used  as  controls. 

Control  (unpulsed)  and  treated  (pulsed)  cells  were  removed 
from  the  electroporation  cuvettes  following  pulsing  and  mixed 
1 : 1  with  a  volume  of  trypan  blue  at  various  times  after 
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pulsing.  A  small  volume  of  this  mixture  was  then  placed  on  a 
hemocytometer,  or  cell-counting  slide.  The  cells  were  then 
counted  as  either  stained  or  unstained.  A  total  of  at  least 
150  cells  were  counted,  and  the  percentage  of  unstained  cells 
was  represented  as  percent  trypan  blue  exclusion. 

We  used  the  trypan  blue  exclusion  assay  as  an  indicator  of 
plasma  membrane  integrity  and  of  cell  viability.  Trypan  blue 
is  normally  impermeant  to  healthy  cells.  When  cell  membrane 
integrity  is  compromised,  the  dye  is  able  to  enter  the  cell  and 
bind  to  protein,  making  the  cell  appear  blue.  Cells  that  take  up 
this  dye  are  usually  considered  dead  or  dying.  However,  in  elec¬ 
troporation  experiments  with  propidium  iodide  as  a  marker  for 
membrane  integrity,  uptake  of  this  marker  after  pulsing  has  also 
been  reported  for  viable  cells  [21].  Because  the  size  of  trypan 
blue  molecules  (960  Da)  is  comparable  to  that  of  propidium 
iodide  (668  Da),  it  cannot  be  assumed  that  observation  of  trypan 
blue  uptake  is  necessarily  an  indication  of  cell  death. 

As  a  more  definitive  means  to  determine  effects  of  sub¬ 
nanosecond  pulses  on  cell  survival,  we  used  a  clonogenic  as¬ 
say.  The  clonogenic,  or  colony-formation,  assay  measures  cell 
survival  based  on  the  ability  of  a  single  cell  to  replicate  success¬ 
fully.  In  this  manner,  it  challenges  every  cell  in  the  population 
to  undergo  continuing  division  [22].  The  clonogenic  assay  mea¬ 
sures  cell  survival  6  days  posttreatment  and  is  linked  to  ultimate 
cell  survival  and  accounts  for  all  lethal  mechanisms  of  cell  in¬ 
jury.  In  this  procedure,  B16  cells  were  placed  in  a  chamber  and 
exposed  to  1000  pulses  at  800  ps  and  550  kV/cm.  After  being 
allowed  to  grow  for  6  days,  single  surviving  clones  had  divided 
to  form  clusters.  Clones  that  had  formed  clusters  of  at  least 
50  cells  were  counted  as  survivors  and  compared  statistically 
to  control  survivors  to  yield  a  survival  fraction  percentage. 

Previous  results  with  60-ns  pulses  indicated  that  cells  died  by 
apoptosis  [9]— [12],  a  natural  mechanism  for  deleting  damaged 
or  unwanted  cells  without  inflammation,  pain,  or  scarring  [23]. 
To  determine  if  apoptosis  mechanisms  accounted  for  cell  death 
after  exposure  to  subnanosecond  pulses,  we  used  a  fluorescent 
cell-permeable  irreversible  inhibitor  of  active  caspases  as  an 
apoptosis  marker  [9]— [12].  Cells  are  incubated  with  the  marker 
and  analyzed  by  flow  cytometry  at  various  times  after  treatment. 
The  flavonoid  phloretin,  which  induced  caspase  activation  by 
a  mitochondria-mediated  apoptosis  mechanism  within  1  h  in 
B16fl0  cells  (data  not  shown),  was  used  as  a  positive  control. 

III.  Results 

The  uptake  of  trypan  blue  was  measured  for  electric  field 
values  between  150  kV/cm  and  1  MV/cm,  depending  on  the 
pulse  number.  Fig.  4  shows  the  results  obtained  for  a  fixed 
electric  field  value  of  550  kV/cm  1  h  after  pulsing.  There  is  no 
trypan  blue  uptake  up  to  about  100  pulses.  With  pulse  numbers 
increasing  beyond  the  threshold,  the  trypan  blue  uptake  (inverse 
of  trypan  blue  exclusion)  increases  exponentially.  This  feature 
is  typical  for  all  the  results  obtained  more  than  1  h  after  pulsing. 

The  trypan  blue  results — obtained  at  1  h  or  more  after 
pulsing — are  consistent  with  the  results  of  clonogenic  assays. 
Results  with  trypan  blue  uptake  and  the  clonogenic  assay  both 
indicated  similar  survival  rates  at  1000  pulses  at  800  ps  and 
550  kV/cm  and  490  kV/cm,  respectively.  This  showed  that  try- 


Fig.  4.  Trypan  blue  exclusion,  1  hour  after  pulsing,  and  caspase  activity, 
measured  4.5  h  after  pulsing,  versus  number  of  800  ps  pulses  with  average 
electric  field  amplitudes  of  550  kV/cm  (trypan  blue  exclusion)  and  490  kV/cm 
(caspase  activity),  respectively. 


Fig.  5.  Temporal  development  of  trypan  blue  exclusion  (625  pulses,  800  ps, 
550  kV/cm). 

pan  blue  uptake  for  these  conditions  can  indeed  be  considered  a 
marker  for  cell  death,  only  if  measured  at  least  1  h  after  pulsing, 
however.  The  mechanism  for  the  cell  death  is,  at  least  in  part, 
caspase-dependent  apoptosis,  which  is  also  shown  in  Fig.  4. 
Cell  survival,  the  inverse  of  cell  death,  decreases  in  the  same 
way  as  the  intensity  of  the  apoptosis  marker  increases. 

The  trypan  blue  exclusion  results  depicted  in  Fig.  4  hold 
for  measurements  where  the  trypan  blue  was  added  to  the 
medium  at  least  1  h  after  pulsing.  An  interesting  result  was 
obtained  when  the  temporal  development  of  trypan  blue  uptake 
was  studied  with  emphasis  on  the  time  before  the  one-hour 
mark.  As  shown  in  Fig.  5,  a  much  larger  percentage  of  the 
cells,  approximately  50%  for  625  pulses  at  550  kV/cm,  take 
up  trypan  blue  right  after  pulsing  (5  min),  compared  to  less 
than  20%,  which  show  trypan  blue  uptake  an  hour  later.  This 
indicates  that  poration  of  the  cell  membrane,  which  is  caused 
by  the  550-kV/cm  pulses,  is  reversible  for  most  of  the  cells 
and  that  the  typical  time  constant  for  pore  closure  under  these 
conditions  is  on  the  order  of  1  h.  The  maximum  appearance  of 
active  caspases,  a  marker  for  apoptosis,  on  the  other  hand,  was 
observed  4.5  h  postpulse. 

For  pulses  with  electric  field  intensity  different  from 
550  kV/cm — varying  between  150  kV/cm  and  1  MV/cm — the 
trypan  blue  exclusion  results  are  similar  to  those  described  for 
the  550  kV/cm  case,  however,  the  number  of  pulses  required 
to  reach  the  threshold  for  trypan  blue  uptake  was  strongly 
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Energy  Density  gE2  tN  (kJ/cm3) 

Fig.  6.  Trypan  blue  exclusion  fraction  versus  electrical  energy  density.  Vari¬ 
ables  were  the  electric  field  intensity  which  was  varied  between  150  kV/cm  and 
1  MV/cm,  and  the  number  of  pulses  which  was  varied  between  1  and  20,000. 

dependent  on  the  applied  electric  field  intensity.  For  150  kV/cm 
pulses  it  took  18,000  pulses  to  reach  a  50%  uptake  mark,  for 
950  kV/cm  pulses,  it  only  required  125. 

The  data  obtained  with  pulses  of  800  ps  duration  seem  to 
indicate  that  trypan  blue  uptake  is  a  dose  effect,  scaling  with  the 
electrical  energy  density.  Fig.  6  shows  the  trypan  blue  exclusion 
fraction  (defined  as:  one  minus  uptake  fraction)  obtained  4  h 
after  N  pulses  were  applied  with  electric  field  amplitudes 
E ,  ranging  from  150  kV/cm  to  1  MV/cm,  plotted  versus  the 
electric  energy  density,  aE2rN ,  where  a  is  the  conductivity 
of  the  suspension  (100  Dcm).  The  exclusion  of  trypan  blue 
decreases  linearly  with  energy  density  above  a  threshold  of 
approximately  0.8  kJ/cm3. 


IV.  Discussion 

A  question  that  we  tried  to  address  with  these  studies  is  about 
the  mechanism  that  causes  cell  death  when  such  ultrashort 
pulses  are  applied.  The  trypan  blue  exclusion  measurements 
indicate  that  the  800-ps  pulses  definitely  have  an  effect  on 
the  outer  membrane.  The  experimental  results  show  immediate 
(minutes  after  pulsing)  uptake  of  the  dye,  which  means  that 
nanopores  with  diameters  large  enough  to  allow  trypan  blue 
(960  Da)  to  pass  have  been  created  in  the  membrane. 

The  effect  on  the  outer  membrane  is  most  likely  a  direct 
effect  caused  by  membrane  charging  and  subsequent  membrane 
poration.  The  notion  of  membrane  charging  followed  by  po- 
ration  of  the  membrane  is  also  supported  by  modeling.  Here, 
this  is  shown  through  molecular  dynamics  (MD)  simulations. 
Details  of  MD  simulations  for  membranes  subjected  to  volt¬ 
age  pulsing  have  been  reported  by  us  elsewhere  [25]— [27] . 
Essentially,  the  MD  technique  mimics  atomistic  many-body 
dynamics  of  a  lipid  bilayer  surrounded  by  hundreds  of  water 
molecules.  Fig.  7  shows  the  simulation  results:  a  2-ns  snapshot 
of  a  section  of  the  cell  membrane  (with  its  molecules  shown 
in  blue)  in  response  to  an  800-ps  560-kV/cm  average  external 
field.  The  time-dependent  intensity  of  the  membrane  field  used 
for  the  MD  simulation  was  obtained  by  applying  the  actual 
external  field  input  as  a  boundary  condition  to  a  distributed 
electrical  representation  of  the  cell  [25],  [26]. 


Fig.  7.  MD  results  for  subnanosecond  pulsing,  with  water  molecules  (shown 
in  red  and  white)  beginning  to  form  nanochannels  through  a  cell  membrane  sec¬ 
tion,  with  its  molecules  shown  in  blue  (snapshot  at  2  ns  after  pulse  application). 

In  Fig.  7,  water  is  represented  by  the  red- white  combination 
(denoting  the  hydrogen-oxygen  atoms)  and  is  present  at  the  top 
and  bottom  of  the  lipid  bilayer.  Small  nanochannels  of  water 
(resembling  “water  nanowires”)  are  clearly  seen.  Evidently, 
a  large  well-defined  channel  is  not  predicted  to  have  formed 
through  the  membrane  due  to  the  ultrashort  duration  of  the 
external  voltage  pulsing.  However,  some  internal  stressing  and 
nanofracturing,  indicative  of  a  shock,  exists  and  portends  the 
eventual  formation  of  aqueous  pores  within  membranes.  The 
salient  feature  is  the  distinct  demonstration  of  electrical  effects 
at  membranes  by  the  subnanosecond  pulse. 

The  modeling  results  on  the  pore  formation  and  the  trypan 
blue  exclusion  studies  (Fig.  5),  which  indicate  pore  closure  with 
a  time  constant  of  about  1  h,  can  be  used  in  a  continuum  model 
to  estimate  the  transport  of  molecules  through  the  membrane 
after  pulsing  (i.e.,  the  rising  part  of  the  trypan  blue  exclusion 
curve  in  Fig.  5).  The  multipulse  model  for  the  pore  trans¬ 
port  model  can  be  constructed  on  the  following  assumptions: 
1)  the  subnanosecond  pulsed  electric  field  is  “on”  for  a  very 
short  time  compared  to  all  other  time  scales;  2)  each  voltage 
pulse  increases  the  size  of  existing  pores  on  the  cell  membrane 
surface  by  an  area,  A S(t);  3)  the  increase  in  pore  area  for  dye 
throughput,  S(t),  after  N  pulses  is  given  as 

s(t)=|]N1/2AS(t)  (1) 

where  the  TV1/2  term  arises  from  the  empirical  scaling  law 
discussed  in  [1].  In  reality,  the  summation  in  (1)  is  an  integral 
over  the  overall  pore  distribution;  and  4)  the  effective  pore 
area  decreases  exponentially  upon  pulse  termination  as  the 
membrane  reseals.  Thus,  the  total,  additional,  pore  area  at  any 
time  t  (assuming  that  the  time  for  the  pulse  application  is  very 
short  compared  to  the  characteristic  time  constant  for  pore 
resealing  r),  is  given  as 

S(t)  =  N1/2ASq  exp(— t/r)  (2) 
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where  A  So  is  the  incremental  area  change  produced  by  a  single 
voltage  pulse. 

The  concentration  of  dye  molecules  inside  cells  increases  due 
to  Fickian  diffusion  through  the  electropores.  Applying  Fick’s 
law,  and  assuming  that  the  dye  concentration  outside  the  cell 
remains  almost  unchanged,  a  time-dependent  equation  for  c(£), 
the  concentration  of  dye  inside  a  cell,  can  be  obtained.  This 
relation,  after  some  algebra,  is 


c(t)  =  Co 


1  —  exp 


(3) 


where  D  is  the  dye  diffusion  coefficient,  d  is  the  membrane 
thickness,  Vc  is  the  average  cell  volume,  and  Co  is  the  initial 
concentration.  In  (3),  T d  represents  a  time  delay  after  voltage 
pulsing  in  introducing  the  dye  into  the  solution.  Also,  time  t  in 
(3)  is  limited  to  durations  during  which  the  resealing  pore  size 
remains  larger  than  the  dye  molecule.  With  trypan  blue  already 
present  in  the  solution  at  the  time  of  voltage  pulsing  of  the  cell 
culture,  Td  =  0.  Integrating  (3),  with  S(t )  given  by  (2),  leads  to 
the  following  relation  for  the  excluded  dye  fraction  /(£,  Td): 


f(t,Td)  =  exp  [{ASoDN^VtVcd)} 

x  {exp(— t/r)  -  exp(-Td/r)} 


(4) 


The  experimental  data  shown  in  Fig.  5,  corresponding  to  an 
average  external  electric  field  of  560  kV/cm  and  625  pulses, 
was  analyzed  based  on  the  above  diffusive  transport  model. 
The  curve  of  Fig.  5  exhibits  saturation  in  trypan  blue  exclusion 
beyond  ^-400  min.  This  implies  that  the  membrane  pores  in 
viable  cells  would  all  have  resealed  to  such  an  extent  that  uptake 
of  the  trypan  blue  molecules  is  no  longer  possible  at  times 
beyond  100  min.  For  such  long  times,  the  fraction  as  a  function 
of  the  Td  can  be  expressed  as 

f(Td)  =  C  exp  [-  {(aSoDN^t)  /(Vcd)}  exp(-Td/r) 

=  Cexp  — K0  exp(— Trf/r )]  (5) 

where 

K0  =  (ASoDN^V)  /(Vc d)  (6) 

and  C  is  a  scaling  constant.  Thus,  there  are  essentially  three 
critical  fitting  parameters,  C,  Kq,  and  r,  for  the  relation  gov¬ 
erning  the  fractional  trypan  blue  exclusion  dependence  on  time 
delay  Td.  For  the  data  given  in  Fig.  5,  the  following  best-fit 
values  were  obtained:  r  =  43.7  min,  K0  =  0.72,  and  C  =  0.87. 

In  the  above,  Ko  is  an  indirect  measure  of  the  relative  ease  for 
trypan  blue  uptake  by  the  cells.  It  is  characterized  by  the  overall 
transport  or  inflow  within  a  given  population  of  cells.  Factors 
that  collectively  influence  its  value  are  pore  area,  diffusion  co¬ 
efficient,  number  of  pulses,  and  membrane  resealing  time.  Also, 
the  Ko  value  can  yield  the  total  pore  area  increment  A  So  per 
pulse  if  the  other  parameters  such  as  the  number  of  pulses,  cell 
volume,  and  dye  diffusion  coefficient  are  known.  In  this  case, 
using  N  =  6 25,  a  cell  radius  of  5  /rm,  r  =  43.7  min  with  D  in 


the  range  1CT11  m2s-1  —  1CT10  m2s-1  yields  ASo  to  be  in  the 
range  of  2.3  nm2-23  nm2 .  The  range  for  D  was  roughly  taken  to 
be  an  order  of  magnitude  lower  than  usually  reported  for  single 
solvated  ions  [28]  due  to  the  large  molecular  size  of  trypan 
blue  and  the  small  openings  of  the  membrane  nanopores.  This 
corresponds  to  radial  increments  in  the  size  of  the  megapore  per 
pulse  in  the  1-3  nm  range,  supposing  that  all  the  ASo  change 
was  to  be  attributed  to  a  dominant  megapore.  In  reality,  there 
will  be  a  distribution  of  pores  with  varying  radii,  and  so  the 
radial  increments  per  pore  and  per  pulse  would  vary  around  a 
value  of  A S/n  where  n  is  the  number  of  pores  in  the  cell. 

Previous  molecular  dynamics  (MD)  calculations  have  shown 
the  pore  radii  in  response  to  nanosecond-pulsed  electric  fields, 
to  be  in  the  1-nm  range  [25],  [29].  Equilibrium  mean  pore  radii 
(prepulse)  are  ro  =  0.4  nm,  and  their  density  is  approximately 
1/fim2  [30],  [31].  In  a  spherical  cell  with  a  diameter  of  10  fim 
we  have,  consequently,  approximately  n  =  1007T  pores.  If  we 
assume  that  each  pulse  contributes  to  an  increase  in  the  pore 
area,  the  mean  (circular)  pore  area  after  N  pulses  would  be 

TIT2  =  TTxl  +  ASoN1/2/!!.  (7) 

For  an  initial  pore  radius,  ro,  of  0.4  nm,  the  pore  number  for  the 
10-/im-diameter  cell  of  n  =  1007 r,  and  ASo  varying  between 
2.3  and  23  nm2  (see  previous  section),  the  radius  r  after  N  = 
625  pulses  is  in  the  range  of  0.47-0.86  nm.  These  values  are 
comparable  to  the  estimated  size  of  nanopores  [25]-[27],  [29]. 

The  results  of  trypan  blue  uptake  studies,  performed  after  ex¬ 
periments  with  various  pulse  parameters,  but  keeping  the  pulse 
duration  constant,  seem  to  indicate  that  the  observed  biological 
effect  is  an  energy  density-dependent  effect  (see  Section  III, 
Results).  However,  this  dependency  on  energy  density  does  not 
hold  true  when  the  pulse  duration  is  varied.  This  is  clearly 
shown  when  we  compare  the  800-ps  data  (Fig.  6)  with  trypan 
blue  uptake  data  obtained  with  pulses  10  ns  long  [Fig.  8(a)]. 
In  Fig.  8(a)  we  have  used  experimental  results  on  trypan  blue 
exclusion  reported  in  [24],  complemented  by  results  of  a  recent 
experimental  study  where  the  procedures  were  kept  exactly 
the  same  as  described  in  [24]  (which  are  identical  to  those 
described  in  this  manuscript).  The  results  show  clearly  that  the 
observed  biological  effect  is  not  an  energy  density-dependent 
effect  but  follows  a  different  scaling  law. 

Such  a  scaling  law  for  biological  effects  of  pulses  with 
durations  short  compared  to  the  plasma  membrane  charging 
time  constant,  but  longer  than  the  dielectric  relaxation  time  of 
the  cytoplasm,  was  discussed  in  [1].  The  law  is  based  on  the 
hypothesis  that  in  this  pulse  duration  range  (which  covers  both 
800  ps  as  well  as  10  ns  pulses),  any  bioelectric  effect  is  due  to 
membrane  charging,  whether  it  is  at  the  plasma  membrane  or 
subcellular  membranes.  It  postulates  that  identical  effects  can 
be  expected  if  the  product  of  applied  electric  field  intensity 
E ,  pulse  duration  r,  and  the  square  root  of  the  number  of 
pulses  N  is  kept  constant.  By  replacing  the  energy  density  in 
the  ordinate  of  Fig.  8(a),  by  this  scaling  parameter,  Er^/N, 
the  trypan  blue  data  for  800  ps  and  10  ns  can  be  described 
by  the  same  parameter,  shown  in  Fig.  8(b).  It  indicates  that  the 
observed  effect  of  trypan  blue  uptake  is  a  primary  effect  caused 
by  membrane  charging  (and  subsequent  poration)  and  is  not  a 
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Fig.  8.  (a)  Trypan  blue  exclusion  fraction  versus  electrical  energy  density  for 

two  pulse  durations:  800  ps  (open  symbols)  and  10  ns  (closed  symbols).  The 
results  show  clearly  that  the  dose  (electrical  energy  density)  required  to  achieve 
identical  effects  (trypan  blue  extension)  is  higher  for  shorter  pulses.  In  this  case, 
for  800  ps  versus  10  ns  pulses,  it  is  approximately  one  order  of  magnitude, 
(b)  Trypan  blue  exclusion  versus  the  scaling  parameter,  Er^/N  for  two 
different  pulse  durations  and  a  wide  range  of  electric  field  amplitudes  and  pulse 
numbers. 


secondary  bioelectric  effect,  e.g.,  caused  by  the  disintegration 
of  the  membrane  after  apoptosis  has  set  in. 


V.  Conclusion 

Whereas  studies  on  nanosecond  pulse  effects  have  already 
led  to  medical  applications,  to  our  knowledge  there  are  no 
experimental  studies  on  the  biological  effects  of  subnanosec¬ 
ond  pulses  with  external  electric  field  amplitudes  greater  than 
10  kV/cm.  By  using  an  800-ps  pulse  generator  we  have,  for  the 
first  time,  focused  on  the  effect  of  these  pulses  on  melanoma 
cells  (B16)  with  trypan  blue  as  an  indicator  of  membrane 
integrity.  We  have  complemented  the  trypan  blue  studies  with  a 
clonogenic  assay  to  determine  cell  death  and  a  caspase  assay  to 
explore  the  mechanism  of  cell  death. 

The  temporal  development  of  the  trypan  blue  uptake  for  625 
and  800  ps,  respectively,  and  550  kV/cm  pulses  shows  that  the 
generated  nanopores  reseal  for  most,  but  not  all,  of  the  cells 
with  a  time  scale  of  40  min  to  1  h.  Approximately  20%  of 
the  cells  are  not  able  to  restore  the  integrity  of  the  membrane 
and  die,  as  shown  by  clonogenic  studies.  This  recovery  of  the 


cell  membrane  is  detrimental  for  applications  such  as  tumor 
treatment,  where  the  goal  is  to  minimize  viability.  However,  for 
applications  where  drugs  are  to  be  delivered  into  cells  through 
the  plasma  membrane,  this  effect  would  be  very  beneficial. 
It  might  allow  delivery  of  certain  drugs  or  nanoparticles  into 
the  cell  without  strongly  affecting  the  viability  of  the  cells. 
In  addition,  the  modeling  studies  indicate  the  possibility  of 
externally  influencing  dye  uptake  in  the  first  hour  after  pulsing 
through  variations  in  the  number  of  pulses  and/or  the  voltage 
magnitude.  This  could  be  important  for  tailored  drug  delivery 
and  selective  gene/molecule  insertion. 

All  our  results  point  to  membrane  charging  and  nanoporation 
[25]— [27],  [29]  as  a  dominant  effect  of  800  ps  pulses,  which  not 
only  causes  reversible  pore  formation  in  the  plasma  membrane 
but  might  also  be  the  reasons  for  irreversible  changes  in  the 
membrane  that  lead  to  cell  death.  That  means  that  the  equivalent 
circuit  shown  in  Fig.  1  is  still  likely  to  describe  the  response  of 
cells  to  electrical  pulses,  rather  than  a  circuit  that  includes  the 
capacitance  of  the  cytoplasm.  This  conclusion  is  supported  by 
results  shown  in  Fig.  8(b),  where  the  trypan  blue  exclusion  is 
plotted  versus  the  scaling  parameter  Er^/N.  Scaling  with  this 
parameter  is  based  on  the  assumption  that  any  observed  effect  is 
caused  by  membrane  charging  [1].  Only  when  this  assumption 
is  valid  do  the  results  of  the  800  ps  experiments  and  the  result 
from  a  10-ns  study  fall  onto  the  same  curve. 

However,  the  fact  that  apoptosis  seems  to  be  the  major 
mechanism  of  cell  death  (Fig.  5)  could  also  indicate  that  direct 
effects  on  membrane  proteins  play  a  role  in  cell  death  and 
that  the  disintegration  of  the  membrane  (observed  1  h  after 
pulsing)  is  partly  a  secondary  effect  due  to  apoptosis.  Because 
we  are  in  a  transition  range  between  membrane  charging  and 
direct  dielectric  effects  with  our  800-ps  pulses,  it  is  likely  that 
both  membrane  charging  and  dielectric  effects  contribute  to  the 
observed  biological  effects.  Further  experiments  with  shorter 
pulses  (250  ps)  are  underway  to  extend  the  parameter  range 
further  toward  the  dielectric  range  and  to  better  discriminate 
between  membrane  charging  and  dielectric  effects,  such  as 
direct  effects  of  the  electric  field  on  membrane  proteins. 

One  reason  to  extend  our  range  of  observation  towards  the 
100-ps  range  is  the  possible  use  of  antennas  as  pulse  delivery 
systems.  The  use  of  needle  or  plate  electrodes  in  therapeutic 
applications  that  rely  on  electroporation  [32]  or  nanosecond- 
pulsed  electric  fields  (nsPEF)  [13]  requires  that  the  electrodes 
are  brought  into  close  contact  with  the  treated  tissue.  This  limits 
the  application  to  treatments  of  tissue  close  to  the  skin,  or 
tissue  close  to  the  surface  of  the  body.  The  use  of  antennas, 
on  the  other  hand,  would  allow  one  to  apply  such  electric  fields 
to  tissues  (tumors)  that  are  not  easily  accessible  with  needles. 
Also,  focusing  electrical  energy  on  the  target  would  reduce  the 
damage  to  the  tissue  layers  surrounding  the  target  and  the  skin. 
The  spatial  resolution  of  an  electric  field  generated  in  tissue 
depends  on  the  pulse  duration  and  the  permittivity  of  the  tissue. 

However,  whereas  the  research  on  wideband  antennas  has 
mainly  focused  on  the  design  and  development  of  antennas 
operating  in  the  far  field,  an  “antenna”  or  better  pulse  delivery 
system  that  is  used  for  bioelectrics  applications  would  need  to 
operate  in  the  near  field.  There  are  several  concepts  for  such 
a  delivery  system.  These  are  mainly  based  on  the  use  of  a 
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Fig.  9.  A  focusing  antenna  using  a  conical  wave  launcher  and  a  prolate 
spheroidal  reflector  [34].  The  electromagnetic  wave  is  launched  at  the  focal 
point  FI,  the  endpoint  of  a  coaxial  cable,  expands  along  the  conical  wave 
launcher,  is  reflected  from  the  prolate  spheroidal  reflector,  and  focused  again 
in  focal  point  F2. 

prolate  spheroidal  reflector,  where  the  pulse  is  launched  from 
one  focal  point  and  reflected  onto  a  second  focal  point.  For 
our  application,  the  second  focal  point  would  be  the  location 
of  a  tumor  inside  the  body.  Much  of  the  design  work  on  such 
delivery  systems  is  done  at  the  Air  Force  Research  Laboratories 
and  at  the  University  of  New  Mexico,  both  in  Albuquerque, 
under  the  guidance  of  C.  Baum  [33]. 

At  Old  Dominion  University,  we  have  concentrated  our 
efforts,  in  cooperation  with  the  team  of  C.  Baum,  on  a  near¬ 
field  antenna  where  the  electric  field  generated  in  the  target  is 
focused  by  means  of  a  prolate  spheroidal  reflector,  in  connec¬ 
tion  with  a  conical  wave  launcher  (Fig.  9)  [34].  The  E-field  at 
the  focal  point  has  only  longitudinal  components.  For  a  100-ps 
risetime,  step-function  input  pulse,  the  focusing  spot  size  (full 
width  at  half  maximum)  in  tissue  is  2-3  cm  along  the  z-axis 
and  <1  cm  along  the  x-axis.  For  a  100-kV  input  pulse,  the  peak 
E-field  at  the  second  focal  point  was  calculated  to  be 
250  kV/cm. 

Our  results,  which  provide  information  on  membrane  in¬ 
tegrity,  and  cell  viability,  indicate  that  the  required  electric 
field  values  for  apoptosis  are  far  above  the  values  that  can  be 
reached  with  wideband  antennas,  even  with  focusing  reflectors. 
However,  by  using  sensitizing  agents  in  combination  with 
wideband  radiation  it  might  be  possible  to  reduce  the  threshold 
for  apoptosis  induction  considerably.  On  the  other  hand,  besides 
apoptosis,  induction  of  nonlethal  effects  such  as  calcium  release 
from  subcellular  structures  [5]-[7]  or  neurotransmitter  release 
can  be  an  attractive  application  of  this  new  method.  Calcium  is  a 
key  regulator  of  numerous  cellular  functions  and  also  influences 
cell  signaling.  This  would,  for  example,  have  implications  for 
electrostimulation.  Because  the  release  of  calcium  has  been 
shown  to  require  much  lower  electric  fields  than  does  apoptosis 
induction,  the  constraints  on  the  pulse  generator  can  be  relaxed 
considerably.  In  general,  the  use  of  subnanosecond  pulses  not 
only  allows  us  to  enter  a  new  field  of  field-cell  interactions  but 
might  open  the  door  to  a  range  of  therapeutic  applications  that 
require  electromagnetic  energy  delivery  into  tissue  not  easily 
accessible  by  solid  electrodes,  such  as  needles. 
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